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Resumen 



El objetivo principal de esta tesis es el desarrollo de un codigo hidrodinamico unidi- 
mensional para modelar curvas de luz bolometricas de supernovas de tipo II plateau 
(SNs II-P) con el proposito final de analizar las propiedades fisicas de los progeni- 
tores de estas supernovas comparando nuestros modelos con una gran base de datos 
de muy buena calidad y cobertura temporal. 

Para poder comparar nuestro modelo con las observaciones hemos derivado calibra- 
ciones para la correccion bolometrica (CB) y la temperatura efectiva basadas en 
fotometn'a BVI . La dispersion tipica de la calibracion encontrada para la CB con- 
duce a una incerteza de 0.05 dex en la luminosidad. 

Usando nuestro codido hemos estudiado en detalle el caso de la SN 1999em, una 
de las SNs II-P mejor observadas, y encontramos un muy buen acuerdo con las ob- 
servaciones usando los siguientes parametros fisicos: E — 1.25 foe, M — 19 M0, 
R — 800 i?0, and = 0.056 M©. En este analisis vemos que es necesario usar 
una extensa mezcla de ^^Ni para reproducir un plateau tan chato como muestran las 
observaciones. 

Hemos calculado curvas de luz bolometricas para todas las SNs de nucstra muestra 
y definido una serie de parametros que nos permiten caracterizarlas. Comparando 
dichos parametros encontramos un rango de 1,15 dex en la luminosidad del plateau, 
un rango dc duracion del plateau de 64 a 103 dias y masas de niquel inferiores a 0, 1 
Mq para todas las SNs en nuestra muestra, con la excepcion de SN 1992am. 
Finamente, hemos calculado una grilla de modelos hidrodinamicos para estudiar la 
dependencia de difcrcntcs obscrvablcs con los parametros fisicos (masas, radio y en- 
crgia dc la explosion) y las corrclacioncs cntrc distintos observables. A partir de este 
analisis vemos que nuestros modelos reproducen la fuerte correlacion observada entre 
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la luminosidad y la velocidad de expansion durante el plateau. 
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Chapter 1 



Introduction 



1.1 Supernova Background 

A supernova ( SN, hereafter) explosion is among the most spectacular events that 
can be observed in the Universe. The term "supernova" is somewhat misleading, 
as such an event does not represent a new star (that is, a "nova"), but instead the 
death of a star. 

The first records on supernovae (SNe) goes back to the second century A.D., 
mainly by Chinese astronomers. They observed and recorded "guest stars" , i.e., 
stars that suddenly appeared in the sky, were visible for a certain length of time, 
and then faded away. The "guest stars" that were visible for a year or longer were 
probably SNe, while the shorter guests were common novae. Among the most fa- 
mous historical SNe are the events of 1054 A.D. (the "Crab Nebula"), 1572 A.D. 
("Tycho's supernova"), and 1604 A.D. ("Kepler's supernova" which appeared only 
few years before the invention of the telescope). Another SN, known as Cas A (Cas 
for Cassiopeia), exploded in our galaxy between 1650 and 1680. Its remnant is a very 
strong radio source but, intriguingly, was not reported by contemporary observers. 
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The occurrence of this phenomenon became more famihar to astronomers at 
the end of the 19th century, with the beginning of the use of photographic plates 
and systematic surveys. Astronomers designated this objects with the genuine name 
"nova". By 1927 ---'lOO novae had been discovered in the magnitude range 16-18. 
However, it was not until the distances to the nearby galaxies were firmly established 
through the observation of Cepheid variables (Hubble, 1925) that the separation of 
ordinary novae and the extraordinary luminous SNe, become clear. 

The first physical explanation for SNe was proposed by Baade and Zwicky 
(1934). They suggested that the SN energy comes from the collapse of ordinary 
stars into neutron stars. This idea is a prescient suggestion for the fate of massive 
stars, but it is not the correct explanation for all types of SNe. Nature has more than 
one way to explode a star. This was revealed clearly by Minkowski (1941) who noted 
that the spectra of SNe indicated at least two different types of objects. One group, 
which he named "Type 1" , was characterized by broad emission features and no signs 
of hydrogen in the spectrum, while the rest, named "Type 11", showed strong, broad 
hydrogen lines dominating the spectrum. Since that time, the classification of SNe 
has evolved parallel to the increasing body of spectroscopy data and our knowledge 
about these objects. In the mid-1980s it was realized that the Type 1 group could 
be further differentiated according to the presence of the Si 11 A6355 line. If this 
line was present in the spectra, the SN was classified as Type la; otherwise it was 
designated as Type lb or Ic, depending on the presence or absence of He lines, 
respectively. Another critical piece of information about the nature of these objects 
was provided by the stellar environments where they were discovered. Supcrnovae of 
types 11, lb and Ic are only observed near star forming regions in late-type galaxies, 
which indicates their association with young stellar environments, while Type la SNe 
(SNe la) occur both in elliptical and spiral galaxies. Currently, it is accepted that 
there are two different physical types: the thermonuclear explosion of a white dwarf 
star in a binary system, and the gravitational collapse of the core of a massive star, 
as suggested by Baade and Zwicky (1934). The latter arc called "core-collapse SNe" 
(CCSNe) and observationally correspond to SNe II, lb and Ic. SNe la correspond to 
the first group of thermonuclear SNe. 

Unlike SNe la, CCSNe form a diverse group in terms of their spectral and 
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photometric properties. SNe II have been subdivided according to their light curve 
shapes into: Il-Plateau (showing plateau-hke hght curves), and II-Linear (showing 
linearly declining light curves). Further subgroups arise from spectroscopic proper- 
ties: Iln (showing narrow emission lines), and lib (showing hydrogen initially but 
afterward evolving to hydrogen- deficient at later times) . There are some SNe Ib/c 
associated with long gamma-ray bursts ( Woosley fc Bloo"ml . 2006 ) which show much 
broader lines than the typical objects of this class. They have been referred to as 
"hypernovae" or broad-lined SNe, due to the large inferred kinetic energi es. Fig- 



Turatto 



ure 11.11 shows a diagram of the current SN classificat ion, as provided by 
(120031 ). A review of SN classification criteria is given by iFilippenkd (119971 ). It should 
be noted that, because of its empirical nature, the current classification scheme is 
prone to evolve as more data are obtained. 
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Figure 1.1: The current classification scheme of supernovae. Type la SNe are associ- 
ated with the thermonuclear explosion of accreting white dwarfs. Other SN types are 
associated with the core collapse of massive stars. CCSNe with explosion energies 
E > lO^^erg are often called hypernovae. 



A typical SN has a luminosity of ^IO^^Lq during a period of weeks to months, 
which is comparable to the luminosity of its host galaxy. This radiative energy 
corresponds to only 0.01% of their typical kinetic energy which is ~10^^ erg (a 
quantity commonly known as 1 foe or, more recently, as 1 Bethe). The kinetic 
energy is essentially determined from the observed expansion velocities, as measured 
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from spectral lines. The characteristic broad P-Cygni profiles in SN spectra, with 
mean expansion velocities of the order to 10^ km s~^, is an indicator of fast-expanding 
atmospheres. In addition, the observed temperatures near maximum brightness are 
~2 Tq, (where Tq stand for the effective temperature of the Sun), which implies 
radii ~10^ R©, assuming spherical symmetry and that the SN radiates as a black 
body. For the case of CCSNe, the kinetic energy represents only 1% of the energy 
released during the explosion. The remaining 99% (~10^^ erg) is carried away by 
neutrinos created during the collapse of the core. Such an explosion mechanism 
was first confirmed with the detection of neutrinos from SN 1987A. This famous 
SN discovered on February 23, 1987 in the Large Magellanic Cloud (LMC) was 
the brightest SN detected since the invention of the telescope and one of the most 
extensively observed stellar objects in the history of astronomy. It was observed in 
all wavelengths from gamma rays to radio wavelength, and for the first time, extra 
solar neutrinos were detected. It presented a peculiar light curve shape which we 
know now that was due to the blue supergiant nature of the progenitor star, later 
identified on photographic plates of the LMC obtained prior to the explosion. 



During the SN explosion a fraction of the progenitor's material is thermonu- 
clearly burned into various elements. Among these products, unstable isotopes of 
iron-group elements play an important role in the subsequent evolution of the object. 
The most abundant radioactive isotope produced is ^^Ni which decays in ^^Co that 
in turn decays to stable ^^Fe. This process provides an additional source of energy 
which is thermalized in the ejecta and later released during a period of weeks to 
months. 



Given their large intrinsic luminosity, SNe have long been considered poten- 
tial probes for extragalactic distance determination and thus for the measurement of 
cosmological parameters that describe the expansion history of the Universe. Histori- 
cally, SNe la have been used for this purpo se due to their high degree of homogeneity. 
Through an empirical correlation found by iPhillipsI ( 1l993l ) between the absolute peak 
luminosity and the width o f the light curve, SNe la can be used to determine distances 
with a precision of ~ 7% ( iHamuy et al.l . Il996[ ). One decade ago, the application of 



this method to high-redshift SNe la in comparison with the low-redshift sample of the 
Calan/Tololo survey led to the remarkable finding of the accelerated expansion of the 
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Universe ( iRiess et all Il998l : iPerlmutter et al.l . Il999[ ). Although, SNe 11-P are fainter 
and less homogeneous than SNe la, they have also been proposed as good distance 
indicators with potential application to cosmology, in an independent way as SNe la. 
The m ost popular methods in the literature are the ex panding photosphere method 



(EPM; 



Eastman et al.l. 



1996 



Dessart fc Hillierl. 



Mitchell et al. 



2nn5a 



ing atmosphere method (SEAM; 
standard- candle method (SCM; iHamuy fc Pinto 



2002 



) , the spectral- fittin g expand- 



Baron et al 



2004h . and the 



2OO2I ) . The former two methods 



require detailed atmospheric modeling, while the SCM is a much simpler approach 
based on an empirical correlation between luminosity and expansion velocity at the 
middle of the plateau phase. In principle, EPM and SEAM provide distances which 
are independent of the standard cosmic "distance ladder" , whereas SCM requires an 
external calibrator in order to provide absolute distances. By comparing low- and 
high-redshift samples, SCM can be potentially employed to determine cosmological 
parameters using data from upcoming SN surveys. 



Recently, a compilation of all SN discoveries within a volume of 28 Mpc during 
a period of tirne of 1 0.5 yr has provided us the relative rates of each sub-type of SNe 
(jSmartt et al.l . 120091 ). This study showed that SNe la represent only 27% of all SNe. 
The rest are CCSNe with the following rates: 59% for SNe II-P, 29% for SNe Ib/c, 
5% for SNe lib, 4% for SNe Iln and 3% for SNe II-L. 



The goal of this thesis is the study of SNe II-P and their progenitors. In § 11.21 
I review our current knowledge on these objects focusing on the motivations for this 



work. In § II. 3[ I summarize the main goals of this thesis and I give a brief overview 
of each of its chapters. 



1.2 Type Il-Plateau SNe 



Type Il-Plateau SNe (SNe II-P) form a well-de fined family characterized by a "plateau" 
in the optical light curve ( iBarbon et al.l . Il979[ l. where the luminosity remains nearly 
constant for a period of ~100 days (see Figure [L2l) . and the presence of prominent 
P-Cygni hydrogen lines in the spectrum (see Figure II. 3p . They are a subclass of 
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the core-collapse SNe (CCSNe) — which includes Type lb, Type Ic, and other sub- 



classes of type II SIN 
greater than 8 Mq (IHeger et al.l 



e — originated by the violent de a th of stars with initial masses 

and sharing, in general 



2003 



Smartt et al. 



I2OO9I ) 



terms, the same explosion mechanism. It has been shown that SNe II-P are the 
most common type of SNe in nature, amounting to 59% of all CCSNe. Apart from 
their astrophysical importance in connection with stellar evolution and the physics 
of the interstellar medium, additional interest on SNe II-P has recently arisen from 
the fact that they have been established as good distance indicators with potential 
application to cosmology, independent of Type la SNe (see § II. ip . 
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Figure 1.2: Light curve SN 2004fx, a typical SN II-P 



Massive stars may suffer considerable mass loss during the early phases of their 
evolution, due to strong stellar winds or transfer to binary companions. Thus, they 
may lose part or all of their outermost envelope of unprocessed hydrogen and helium. 
The vast diversity observed among CCSNe is related with the properties of the 
progenitor star. In the current picture, SNe II result from progenitors which are 
able to retain a significant fraction of their external hydrogen layers and possible 
have the least massive progenitor of all CCSNe. This picture is combined with 
hydrodynamical models of SNe II-P, which show that a red supergiant progenitor 
with an exte nsive H envelope is necessary in order t o reproduce the platea u-shaped 



light curves (iGrassberg et al 



1971 



Falk fc Arnett 



1977; 



Chevalierl . Il976l ) . Recent 
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Figure 1.3: Spectra of SN 2004fx, a typical SN II-P 



direct detection of the progenitors of several SNe II-P sugg est in fact that SNe II - 
P arise from red s uperg iant stars with Mz ams < ISMq (I Van Dyk et al.l (120031 ): 



Smartt et al. 



mm and ISmartt et all mm ). 



There is general agreement that the explosion of a massive star is originated by 
the collapse of its central parts into a neutron star or black hole when the iron core 
is formed at the end of the star evolution, and further nuclear burning no longer pro- 
vides thermal pressure to support the star. During the explosion, the heavy elements 
(mostly a-nuclei) synthesized inside the SN progenitor over its life are spread at very 
high speed into the interstell ar medium. Despite t he intensive theor etical modeling 
done during recent years (see 



Burrows et al. 



2006 



Janka et al. 



2007. and references 



therein) the mechanism that transfers the energy released during the collapse of the 
core to the envelope, still remains unknown. The approach usually adopted to model 
CCSNe is to decouple the explosion into two independent parts: a) the core collapse 
and the subsequent formation of a shock wave, and b) the ejection of the envelope. 
Based on the analysis of the propagation of the shock wave through the envelope, 
independently of how the shock is formed, it is possible to study the observational 



outcome of the e xplosion such as light curves and spec tra 



extensively used (IFalk &: Arnett 



1977 



Grassberg et al 



1971 



lis approach h as been 



Woosley 



19881 . among 



others) and has led to the conclusion that the main factors influencing the outcome 
are the explosion energy and the progenitor properties (mass and radius). 
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Observationally, SNe II-P show a wide range of plateau lumin osities (Lp) and 
durations (At„), exp ansion velocities (vexp), and nickel masses (Mnj) ( jYoung fc Branch! . 



1989 



Hamuyl . l200l[ ). The morphology of the light curve (LC) of SNe II-P has been 



shown to be connected with physical properties of the progenitor object such as 
ejected mass (M), explosion energy (E), and pre-supernova radius (R). The relations 
between these physical and observable parameters were first derived analytically by 



Arnett 



(Il980[ ) and then gener alized bvlPopov 



19931) . Numerical calibrations of these 



relations were then given by iLitvinova fc NadezhinI (119831 119851 ) (LN83 and LN85 
hereafter, respec tively) based o n a g rid of hydrqdynan iical models for different values 



of M, R and E. 



Hamm 



vl (I2OO3I ) 



and 



NadyozhinI (120031 ) applied such calibrations to a 



set of ~20 SN II-P and thereby derived masses, radii and explosion energies for their 
sample. However, these studies have not been fully satisfactory owing to 1) the lack 
of good-quality data, 2) the use of simplified relations between ill-defined and hard- 
to-measure photometric and spectroscopic parameters, and 3) the fact that some of 
the models are based on simplified physical assumptions. Some of the weaknesses of 
the LN83 and LN85 models are that they did not include the effect of nickel heating 
in their calculations, the use of old opacity tables, the neglect of any effect from line 
opacities, and the simplified initial pre-supernova models adopted. 



This thesis was motivated by the need to remedy such problems in order to 
improve our knowledge of the physical properties of SN II-P progenitors. Further- 
more, the avail ability of an enlarge d data set of BVI light curves and spectra of 
~30 SNe II-P (IHamuy et al.l . I2OIOI ) justifies a detailed study. To accomplish our 
purpose, we have developed our own hydrodynamical model (IBersten et al.l . l2010h 
which allows us to perform a comparison between models and data in a consistent 
way. 



During the develo pment o: 



:his t hesis there have been important advances in 
this field. For example, lUtrobinI ( 120071 ) performed a detailed analysis of SN 1999em 
and also studied how several physical parameters affect the LC. In that work, the 
author provided relations between ph ysical and observed param eters for SNe II-P 
similar to SN 1999em. More recently, iKasen fc Woosleyl (|2009[ ) calculated a set of 
model LC and spectra of SNe II-P, for different masses, metallicities, and explosion 
energies. They used their models to describe the dependence of plateau luminosity 
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and duration on explosion energy and progenitor mass. Nevertheless, the relations 
they found are simple and easy to apply only in the extreme case of no ^^Ni pro- 
duction. When ^^Ni is considered, the relations involve more parameters, which 
hinders their applicability to obtain physical parameters from observations. Other 
works have also explored the effect of sever al physical paramet e rs on the LC and 



Youn 



20041 ) ■ al- 



other observational properties of SNe II-P (iChieffi et al.l . 120031 : 
though without attempting to derive general relations. In recent years, there has 
been a number of works which analy zed the physical properties of individual SNe II 
P based on hydrodynamical models (IBaklanov et al 
among others). 



2005 



Utrobin fc Chugai 



2009 



In addition to the comparison with hydrodynamical models, there is an alter- 
native way to derive progenitor masses of SNe II-P, namely, using pre-supernova 
images to search for the progenitor star and measuring its brightness and color to 
derive a mass in connection with a stellar evolution model. The last decade has wit- 
nessed direct discoveries of several SN progenitors providing a remarkable and rapid 
progress in this field. Currently, there are three SNe II-P (SN 2003gd, SN 2005cs 
and SN 2008bk) with a firm detection of th e progenitor and ~17 others without 
a positive detection (see 



Smartt et al. 



(120091 ) for a review). The latter cases are 
still useful since they provide upper limits to the progenitor masses. These studies 
have suggested a minimum of the stellar mass for SNe II-P of 8 ± IM© and a max- 
imum of 16.5 ± 1.5Mq. Although there seems to be consensus on the lower mass 
limit of mass for SNe II-P, the upper limit derived by this method is systematically 
lower than the pro g enitor mas ses deriy e d by modelling the light curves, as noted by 



(2008 


) and 


Smartt 


(2009) 



P with hydrodynamical masses which have been studied in pre-supernova imaging 
(namely, SN 1999em, SN 2004et and SN 2005cs) and a few more that have been 
studied using semi-analytical models. In all cases the mass estimated by hydrody- 
namical models is higher than the estimate or upper limit given by pre-supernova 
imaging. This discrepancy poses a very interesting and unsolved problem which is 
necessary to address. 
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1.3 This thesis 



Supernovae are very relevant astrophysical objects: they affect the energetic and 
chemical evolution of galaxies, and are also important distance indicators for cos- 
mology. Determining what ranges of physical parameters (masses, radii and energies) 
give rise to each type of CCSNe is still an open question. One way to estimate such 
parameters is by comparing observations (light curves, colors, spectra) with hydro- 
dynamical models. In addition to this, the LC provides a useful tool to test models of 
stellar evolution of massive stars and unveil the internal structure (chemical compo- 
sition, mixing and mass distribution) of the progenitor object. Furthermore, models 
may provide a theoretical explanation for the correlations that are employed in the 
determination of distances to SNe II-P, and guide future refinements of this tech- 
nique. 

The main goal of this thesis is the development of a one-dimensional, Lagrangian 
hydrodynamic code to model bolometric light curves of SNe II-P with the ultimate 
purpose of analyzing the physical properties of an available large sample of SNe with 
highly-precise, well-sampled observations. 

I have organized the thesis as follows. In Chapter [21 I give a description of 
the sample of ~30 SNe II-P used in this study, the sources of the data, corrections 
for extinction and distances. In Chapter [3l I present our one- dimensional, flux- 
limited, Lagrangian hydrodynamical code for the modeling of bolometric light curves 
of SNe II-P. The micro-physics and pre-supernova models used in our calculation are 
also described in Chapter [3l More technical details of the calculations are given in 
appendices |A] and [HI In Chapter [H I explain the method used to compute bolometric 
corrections (BC) for SNe II-P which are necessary to derive bolometric light curves 
from BVI photometry, this allows us to compare our model with observations. In 
Chapter [5l I analyze the case of the prototype SN 1999em, and discuss how our 
model compares with previous hydrodynamical studies of this object. In Chapter [6l 
I derive bolometric LC for our sample of SNe and derive physical parameters for 
this sample. I also study the correlations between different observables. Finally, in 
Chapter [TJ I summarize the main conclusions of this work. 
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Chapter 2 

Sample of supernovae 



One of the main motivations for this thesis is the availabihty of an unprecedented 
data set of BVI hght curves and spectra of ~35 SNe which offers the opportunity to 
significantly improve our understanding of SNe II-P and their progenitors. This sam- 
ple of SNe II-P was observed in the course of four systematic follow-up programs: the 
Cerro Tololo Supernova Program (1986-2003), the Caldn/Tololo Supernova Program 
(CT; 1990-1993), the Optical and Infrared Supernova Survey (SOIRS; 1999-2000), 
and the Carnegie Type II Supernova Program (C ATS; 2002-2003). Ad ditionally, I 
i ncluded four SNe from the liter ature: SN 1999gi Jheonard et al.l . 120021). SN 2004di 



Jvinko et al.lboOGk SN 200 4et jsahu et al.l . bo06l ). and SN 2005cs jPastorello et al. 



200' 



by 



Tsvetkov et 



Van Dyk et al 



al. 



(2o3 



20061). Complementary pho tometry for SN 2003gd obtained 



and 



Hendry et al.l (|2005[ ) was also incorporated. A com- 



plete listing of the SNe used in this work is given in Table 12.11 where I present the 
SN, host galaxy name, the SN equatorial coordinates, the heliocentric redshift of the 



host galaxy and their sources, the extinction due to our own Galaxy ( ISchlegel et al 



19981 ) and the corresponding reference for the photometric and spectroscopic data of 



each supernova. 
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Table 2.1: List of Type II Supernovae 



C T\T 1\T ^ 1^ ^ 


Host Galaxy 


a 


(J2000) 


o(J2000) 






(s) 


E{B - 




Hererences 


1991al 


LEDA 140858 


19 


42 


24, 


.00 


-55 06 23.0 





.01525 


HP02 


0. 


051 


1 


1992af 


ESO 340-G038 


20 


30 


40, 


.20 


-42 18 35.0 





.01847 


NED 


0. 


052 


1 


1992am 


MCG -01-04-039 


01 


25 


02, 


.70 


-04 39 01.0 





.04773 


NED 


0. 


049 


1 


1992ba 


NGC 2082 


05 


41 


47, 


.10 


-64 18 01.0 





.00395 


NED 


0. 


058 


1 


1993A 


anonymous 


07 


39 


17, 


.30 


-62 03 14.0 





.02800 


NED 


0. 


173 


1 


1999br 


NGC 4900 


13 


00 


41, 


.80 


4-02 29 46.0 





.00320 


NED 


0. 


024 


2 


1999ca 


NGC 3120 


10 


05 


22, 


.90 


-34 12 41.0 


0, 


.00931 


NED 


0. 


109 


2 


1999cr 


ESO 576-G034 


13 


20 


18, 


.30 


-20 08 50.0 





.02020 


NED 


0. 


098 


2 


1999em 


NGC 1637 


04 


41 


27, 


.04 


-02 51 45.2 





.00267 


NED 


0. 


040 


2 


1999gi 


NGC 3184 


10 


18 


17, 


.00 


-1-41 25 28.0 


0, 


.00198 


NED 


0. 


,017 


3 


200210 


MCG -1-00-03-054 


01 


01 


16, 


.80 


-01 05 52.0 





.05140 


NED 


0. 


036 


4 


2002fa 


NEAT J205221.51 


20 


52 


21, 


.80 


-1-02 08 42.0 





.06000 


NED 


0. 


099 


4 


2002gw 


NGC 922 


02 


25 


02, 


.97 


-24 47 50.6 





.01028 


NED 


0. 


020 


4 


2002hj 


NPMlG-l-04.0097 


02 


58 


09, 


.30 


-1-04 41 04.0 





.02360 


NED 


0. 


115 


4 


2002hx 


PGC 23727 


08 


27 


39, 


.43 


-14 47 15.7 





.03099 


NED 


0. 


054 


4 


2003B 


NGC 1097 


02 


46 


13, 


.78 


-30 13 45.1 





.00424 


NED 


0. 


027 


4 


2003E 


MCG -4-12-004 


04 


39 


10, 


.88 


-24 10 36.5 





.01490 


J09 


0. 


048 


4 


ZUUO 1 




no 


14 


11, 


.06 


-LI AA zis n 


n 
u 


09701 




n 
u. 


(Til 


A 


2003bl 


NGC 5374 


13 


57 


30, 


.65 


-1-06 05 36.4 


0, 


.01459 


J09 


0. 


,027 


4 


2003bn 


2MASX J10023529 


10 


02 


35, 


.51 


-21 10 54.5 


0, 


.01277 


NED 


0. 


,065 


4 


2003ci 


UGC 6212 


11 


10 


23, 


.83 


-1-04 49 35.9 


0, 


.03037 


NED 


0. 


,060 


4 


2003CI1 


IC 849 


13 


07 


37, 


.05 


-00 56 49.9 


0, 


.01811 


J09 


0. 


,021 


4 


2003CX 


NEAT J135706.53 


13 


57 


06, 


.46 


-17 02 22.6 


0, 


.03700 


NED 


0. 


,094 


4 


2003ef 


NGC 4708 


12 


49 


42, 


.25 


-11 05 29.5 


0, 


.01480 


J09 


0. 


046 


4 


2003fb 


UGC 11522 


20 


11 


50, 


.33 


-1-05 45 37.6 


0, 


.01754 


J09 


0. 


,183 


4 


2003gd 


M74 


01 


36 


42, 


.65 


-1-15 44 20.9 


0, 


.00219 


NED 


0. 


,069 


4,5,6 


2003hd 


MCG -04-05-010 


01 


49 


46, 


.31 


-21 54 37.8 


0, 


.03950 


NED 


0. 


,013 


4 


2003hg 


NGC 7771 


23 


51 


24, 


.13 


-1-20 06 38.3 


0, 


.01427 


NED 


0. 


,074 


4 


2003hk 


NGC 1085 


02 


46 


25, 


.76 


-1-03 36 32.2 


0, 


.02265 


NED 


0. 


,037 


4 


2003hl 


NGC 772 


01 


59 


21, 


.28 


-1-19 00 14.5 


0, 


.00825 


NED 


0. 


,073 


4 


2003hn 


NGC 1448 


03 


44 


36, 


.10 


-44 37 49.0 


0, 


.00390 


NED 


0. 


,014 


4 


2003ho 


ESO 235-G58 


21 


06 


30, 


.56 


-48 07 29.9 


0, 


.01438 


NED 


0. 


039 


4 


2003ip 


UGC 327 


00 


33 


15, 


.40 


-1-07 54 18.0 


0, 


.01801 


NED 


0. 


066 


4 


2003iq 


NGC 772 


01 


59 


19, 


.96 


-1-18 59 42.1 


0, 


.00825 


NED 


0. 


,073 


4 


2004dj 


NGC 2403 


07 


37 


17, 


.00 


-1-65 35 58.1 


0, 


.00044 


NED 


0. 


,040 


7 


2004et 


NGC 6946 


20 


35 


25, 


.30 


-1-60 07 18.0 


0, 


.00016 


NED 


0. 


,342 


8 


2005CS 


NGC 5194 


13 


29 


53, 


.40 


-1-47 10 28.0 


0, 


.00154 


NED 


0. 


,035 


9,10 



^ Heliocentric host-galaxy redshifts 

Sources of host-galaxy redshifts. HP02: iHamuv fc Pintol ll2002t) : J09: Ijones et all l l2009l) : NED: 
NASA/IPAC Extragalactic Database 

Ref erences. (1) C a ldn/T ololo Supernova Pro gram; (2) S OIRS; (3) [Leonard et al.l ll2002fl: (4) CATS ; 
(5 ) IVan Dvk et al.l ll2003f) : (6 ) iHendrv et al.| lf2005i) : (7) IVinko et al.l l|2006l ): (g) ISahu et al.l l|2006i) : 
f9 lPastorello et al.l ll2006t) : f 10) iTsvetkov et al.l ll2006h 
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2.1 Observations 



The observations were made with telescopes from Cerro Tololo Inter- American Ob- 
servatory (CTIO), Las Campanas Observatory (LCO), the European Southern Ob- 
servatory (ESO) in La Silla and the Steward Observatory (SO). Table 12.21 shows 
the telescopes and instruments used to obtain the photometry and spectroscopy of 
this data set. CCD d e tectors and standard J o hnson -Kron-Cousins UBVRIZ filters 
(iJohnson et al.l . llQGGt ICousinsl . llOTlt iHamuy . were employed in all the cases, 



and for a small subset of SNe, observations in the JHK filters were also obtained. 
Additionally, low resolution optical spectra were taken for each SN at various epochs 
using CCD detectors in combination with different gratings/grisms and blocking fil- 
ters. Currently, all of th e optical data and sp ectra have been reduced and they are 
in course of publication (IHamuy et al.l . I2OIOI ). The data reductions were performed 
using IRAFlj procedu res including, among others, the host galaxy subtraction (for 
more details, see also iHamuyl . I2OOII ). 



2.2 AKA Corrections and Distances 



In order to study the physical properties of our SNe and to compare them with each 
other it is necessary to correct the observed magnitudes by host-galaxy extinction 
(A/iost), redshift (K-correctio n), and Galacti c extin ction (Ac) dubbed AKA correc- 
tions as described in detail by 



Olivares et al. 



( I2OO9I ). It is also important to place all 



of the SNe in the same distance scale by correcting for systematic differences among 
the different distance estimate methods. 



The AKA corrections and distances for most SNe were calculated by 



Olivares et al, 



(120091 ). While th e determination of Ga lactic extinction is straightforward using the 
IR dust maps of ISchlegel et al.l (1l998l l and assuming a standard reddening law as 



^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by 
Association of Universities for Research in Astronomy, Inc., under cooperative agreement with the 
the National Science Foundations 
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Table 2.2: Telescopes and instruments used in the photometric and spectroscopic 
observations 



Telescope 


Instrument 


Spec/Phot'' 


CTIO 0.9m 


CCD 


P 


YALO 1.0m 


ANDICAM 


P 


YALO 1.0m 


2DF 


S 


CTIO 1.5m 


CCD 


P 


CTIO 1.5m 


CSPEC 


s 


Blanco 4.0m 


CSPEC 


s 


Blanco 4.0m 


2DF 


s 


Blanco 4.0m 


CCD 


p 


Swopc 1.0m 


CCD 


p 


du Pont 2.5m 


WFCCD 


S/P 


du Pont 2.5m 


MODSPEC 


s 


du Pont 2.5m 


2DF 


s 


du Pont 2.5m 


CCD 


p 


Baade 6.5m 


LDSS2 


S/P 


Baade 6.5m 


B&C 


s 


Clay 6.5m 


LDSS2 


S/P 


ESO 1.52m 


IDS 


s 


Danish 1.54m 


DFOSC 


S/P 


ESO 2.2m 


EFOSC2 


s 


NTT 3.58m 


EMMI 


s 


ESO 3.6m 


EFOSC 


s 


Kuiper 61" 


CCD 


p 


Bok 90" 


B&C 


s 



Whether the instrument was used for photometry (P), spectroscopy (S) or both (S/P). 



14 



3.1, determin i ng th e extinction due 



Olivares et al. 



fcoOOh addressed this 



that given by ICardeUi et al.l (Il989[ ) , with Ry 
to host-galaxy dust is a more difficult task, 
issue by assuming that all SNe II-P should evolve to a similar asymptotic intrinsic 
color at the end of the plateau phase. This is expected because of the hydrogen 
recombination nature of the photosphere. Therefore, color indices measured at the 
end of the plateau can be used to determine reddening by comparing with a sample 
of unreddened SNe II-P. 



The distances to this sample are determined using the Standardized Candle 
Method (SCM), with the exception of SN 2003ef and SN 2005cs. For SN 2003ef I 
adopte d the EPM distance es timated by 



Jones et al 



scale of lOlivares et al 



(120091), converted to the distance 



( 120091 ) using the conversion coefficients given by the authors. 
The resulting value is 48.55 ± 18.4 Mpc. For SN 2005cs, I used the distance modulus 
of = 29.62 ± 0.25 mag given by Pastorello et al 2006, which corresponds to 8.4 ± 1 
Mpc. Table 12.31 shows host-galaxy extinctions and distance moduli for our sample 
of SNe. 



2.3 Time Reference Frame 



An additional important step which is necessary in order to compare the data of 
different SNe is to place them in a consistent reference frame of time. This is achieved 
by establishing an origin of time for each SN which places time intervals for the 
observations of all SNe on uniform basis. Time intervals are then corrected for time 
dilation using the redshift of each SN. 



A natural choice of the origin of time would be the moment of the explo- 
sion. However, there are only 12 SNe in the sample with an accur ate estimate 



of th e explosion time (to) based on EPM and pre-explosion imaging ( jjones et al. 



2009|). These estimates were obtained using two inde pendent atmospheric models 



by 



Eastman et al. 



(119961) and iDessart fc Hilhed fl2005af ) (E96 and D05 hereafter, re- 
spectively), which led to two different values of to for each SN. In Table l273l I show 
both values of to for those 12 SNe II-P. I use additional information of the explo- 
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Vinko et al. 



Smartt et al. 



sion time found in the literature for SN 1 991al, SN 199 af, SN 1992am, SN 1993A 
SN 1999ca, and SN 1999cr (see Table 2 of 
to = 471 6.5 ±21 given b^^ 
given by 



Hamuvl . [20031). For SN 2003gd a value of 
(120041 ) ■ for SN 2004dj a value of to = 5167±21 



mm . for SN 2004et a value of to = 5270.5 : 
2006[) and finally for SN 2005cs a value of to = 5549.9 ± 1 given by 
(j2006|). All these values of to are given respect to JD = 2448000. 



2 fISahu et al. 



Pastorello et al. 



An alternative approach was introduced by lOlivares et al.l ( l2009l ). In such work, 
the origin of time was defined as the moment at the middle point of the transition 
between the plateau and the radioactive tail (tp^)- The values of tpx for the sample 
of SNe are given in Table 12.31 With this definition, I was able to place all the SNe 
of our sample in a common time frame. Both origins of time, to and tpT, are used 
in this thesis, depending on the context of the analysis. 
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Table 2.3: Host-galaxy extinctions, distances and explosion times for the sample of 
SNe II-P 











40(005)" 




SN Name 






[JD-2448000] 


[JD-2448000] 


[JD-2448000] 


1991al 


-0.17(21) 


34.05(28) 






521.998 


1992af 


-0.37(21) 


34.15(28) 






858.583 


1992am 


0.52(23) 


36.57(025)<* 






947.43 


1992ba 


0.30(21) 


31.33(28) 


883.9(3) 


879.8(5.6) 


1013.714 


1993A 


0.06(25) 


35.59(04)<i 






1106.96 


1999br 


0.94(25) 


31.86(31) 




3275.6(7.7) 


3395.508 


1999ca 


0.25(21) 


33.29(11)<^ 






3373.19 


1999cr 


0.12(21) 


34.56(28) 






3348.88 


1999em 


0.24(21) 


30.00(28) 


3476.3(1.1) 


3474.0(2) 


3590.136 


1999gi 


1.02(21) 


30.50(28) 


3517.0(1.2) 


3515.6(2.4) 


3647.759 


0210 


0.31(36) 


37.31(31) 






4591.917 


2002fa 


-0.35(26) 


37.16(29) 






4576.632 


2002gw 


0.18(22) 


33.41(28) 


4557.9(2.7) 


4551.7(7.6) 


4662.195 


2002hj 


0.24(22) 


34.93(28) 






4660.231 


2002hx 


0.38(22) 


35.49(28) 






4657.707 


2003B 


-0.09(21) 


32.18(27) 






4713.878 


2003E 


0.78(23) 


34.01(28) 






4765.663 


2003T 


0.35(21) 


35.17(28) 


4654.2(2.7) 


4648.9(3.4) 


4759.381 


2003bl 


0.26(21) 


34.07(30) 




4692.6(2.8) 


4804.609 


2003bn 


-0.04(21) 


33.60(28) 


4693.4(2.7) 


4687.0(9) 


4817.0 


2003ci 


0.78(23) 


35.30(38) 






4817.720 


2003cn 


-0.04(23) 


34.81(28) 






4804.986 


2003CX 


-0.27(25) 


36.20(29) 






4830.927 


2003ef 


0.98(21) 


33.43(1.2)'= 


4759.8(4.7) 


748.4(15.6) 


4869.095 


2003fb 


1.24(23) 


34.49(29) 






4874.252 


2003gd 


0.33(21) 


29.98(28) 






4840.108 


2003hd 


0.01(21) 


35.86(28) 






4953.053 


2003hg 


1.97(24) 


33.31(28) 






4998.604 


2003hk 


0.44(25) 


34.41(28) 






4957.908 


2003hl 


1.72(23) 


32.04(28) 


4872.3(1.7) 


4865.4(5.9) 


5006.708 


2003hn 


0.46(21) 


31.13(27) 


4859.5(3.8) 


4853.8(9.3) 


4962.898 


2003ho 


2.19(21) 


34.13(28) 






4919.688 


2003ip 


0.56(22) 


33.76(28) 






5002.872 


2003iq 


0.25(22) 


32.43(28) 


4909.6(4.3) 


4905.6(9.5) 


5019.625 


2004dj 


-0.09(23) 


28.14(29) 






5286.428 


2004et 


0.13(27) 


28.37(31) 






5389.0 


2005CS 


0.72(30) 


29.62(0.25)'' 






5665.705 



^ Host-galaxy extinctions and distance moduli calculated bv lOlivares et aTl ||2009|) . 

^ Explosion times calculated by I Jones et al.l ||2009| ) using two atmospheric models (E96 and D05). 

*= Times at the middle point of the transition between the plateau and the radioactive tail, as calculated by 
lOlivares et"al] | |2009|) . 

Redshift adopting = 68. 

EPM distance estimated by I Jones et al.l l|2009h . converted to the distance scale of lOlivares et al" ] l|2009r i. 
' Distance modulus given by Pastorello et al 2006. 
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Chapter 3 



Hydrodynamical Code 



The main part of my thesis work is devoted to develop a one-dimensional hydrody- 
namic code in order to model bolometric light curves of SNe II-P. Due to the difficulty 
of modelling the core collapse from first principles, the approach extensively used in 
the literature, and followed in this work, is to decouple the SN explosion problem in 
two independent parts: a) the core collapse and formation of the shock wave (SW), 
and b) the ejection of the envelope. Although the exact nature of the mechanism 
that rips the object in two parts is not clear, the great differences in energetics and 
time scales (seconds for the collapse and days for the ejected envelope) allow the the- 
oretical description of a SN II explosion to be separated into internal and external 
problems. The energy transferred to the envelope (which we called explosion en- 
ergy) plays the role of a coupling parameter between the internal and external parts 
of problem. In addition, the processes which control the envelope ejection and the 
supernova radiation do not depend on how the energy is transferred to the envelope 
as long as this process occurs in a short enough time. Based on the propagation 
of the SW through the envelope, independently of how the shock is formed, it is 
possible to study the observational outcome of the explosion, such as light curves 
and spectra. 

In § 13. H I describe the numerical method, equations, micro-physics and ini- 
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tial models used to compute bolometric light curves of SNe II-P. That Section is 
complemented by Appendices |A] and [B] which provide a technical description of the 
methods. The structure of the code is presented in § 13.21 A discussion of the adopted 
approximations is given in § 13.31 



3.1 Calculation Method 



Our supernova models are computed by numerical integration of the hydrodynamical 
equations assuming spherical symmetry for a self-gravitating gas. Radiation trans- 
port is treated in the dif f usion approximation with the flux-limited prescription of 



Levermore fc PomraningI (119811 ). The explosion is simulated by injecting a certain 
amount of energy near the center of the progenitor object during a very short time as 
compared with the hydrodynamic time-scale. This energy induces the formation of 
a powerful shock wave that propagates through the progenitor transforming thermal 
and kinetic energy of the matter into energy that can be radiated from the stellar sur- 
face. To calculate shock waves, we include, as is usually done, an artificial-viscosity 
term in the equations of moment and energy. 

The equations and numerical method used are discussed in § 13.1.11 In § 13.1.21 
I give a description of the constitutive relations included in the code. The energy 
deposited by radioactive decay is discussed in § 13.1.31 Finally, in § 13.1.41 the initial 
models are described. 



3.1.1 Equations for radiation transport and hydrodynamics 

The differential equations that describe the hydrodynamics and radiative transfer of 
the system assuming spherical symmetry and diffusion approximation in Lagrangian 
coordinates are: (a) velocity definition, 
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(b) mass conservation, 



(c) momentum conservation, 



dr 

^ = «. (3.1) 



V = i = -?^, (3.2) 



p 3 dm ' 



(d) energy conservation. 



| = + (3.3) 

at am 



and (e) radiative energy transport. 



^ = -(4-^)^ ^ f^, (3.5) 
3k am 



where m is the Lagrangian mass coordinate chosen as independent variable in place 
of the radius (r) to describe the structure of the object, u is the velocity, V is the 
specific volume, P is the total pressure (of gas and radiation), and q is the artificial 
viscosity which is included in the equations to spread the pressure and energy over 
several mass zones at the shock front. There are many expressions for the artificial 
viscosity, all dependent on the velocity gradient, which aim at providing a convenient 
interpol ation scheme between unshocked an d shocked fluid. We adopt the expression 



given by I Von Neumann fc Richtmyerl ( 119501 ) (see Appendix E is the total internal 



energy per unit of mass, including gas and radiation, eNi is the energy deposited by 
the radioactive decay of nickel as I describe in § 13.1.31 In equation (13. 4 p we do 
not consider other sources of cooling or heating, such as losses due to neutrino 
processes or energy released by thermonuclear reactions, because the energy of the 
SW overwhelms these other sources. Even if neutrinos are very important in the 
formation of the shock wave as the explosion depends noticeably on the efficiency 
of their energy transfer, most of them are emitted before the shock wave reaches 
the stellar photosphere, so they have no effect on later epochs of the SN evolution 
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see 



Hillebrandt . 


1990; 


Burrows, 


1991; 


Janka et al.. 


2007) 



20071 ). The energy released by 



explosive nucleosynthes is is much less than the ener gy of t he shock wave as has been 
previ o usly est ablis hed 



Imshenik fc Nadezhin (see 



19651); IWooslev fc Weaver 



^ (see 

1990 ); ArnettI (see 1996 )]. T is the temperature of both matter and radiation, n is 
the Rosseland-mean opacity and L is the luminosity. Finally, A is the so-called "flux- 
limiter", included in the diffusion equation to ensure a smooth transition between 
diffusion and free-streaming regimes to assure causality. The expression adopted for 
A is. 



6 + 3i? 



where 



Note that when — )■ oo (short mean free path), i? — )■ and A — ?■ 1 as expected for 
the diffusion limit. But as — )■ (long mean free path), i? — )■ oo and A — )■ 3/-R. 
Thus, it can be seen that the flux is limited to the value caT"^, as required physically. 
However, fiux-limited diffusion is not the correct solution for radiation transfer in 
optically thin regions; it is simply an interpolation between the physically correct 
optically thick and thin limits. 

The quantities P, and k, are functions of p, T, and chemical composition. 
Further details on these constitutive relations are given in § 13.1.21 

The overall problem is to solve the complete structure of the object, i.e., de- 
termining the dependent variables r, u, V, T, and L as a function of m and time. 
Therefore, the problem reduces to solving a system of five coupled partial differential 
equations (I3.ip - (I3.5I) together with the constitutive relations (see § I3.1.2p and the 
expression for q, A and eNi (see § I3.1.3p . In addition, it is also necessary to provide 
boundary and initial conditions. The boundary conditions that we have adopted are: 
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u = near the center (at m = Mcore, where we typically set Mcore = 1-4 M©), and 
-Pgas = p = at the surface (m = M). As initial condition, we have used different 
initial models (see § I3.1.4p . in hydrostatic equilibrium. As discussed above, we sim- 
ulate the explosion by artificially adding internal energy almost instantaneously in 
the central region of the core. 



We have developed a code in FORTRAN to numerically solve this system of 
differential equations using the method of finite differences. As a first step, it is 
necessary to discretize each equation in order to transform the system of differen- 
tial equations into an algebraic one. It is important to note that the discretization 
process is not trivial as t here is no unique way to discretize a differential equation 
(iBowers &: Wilsonl . Il99l[ ). Details of the discretization and the scheme of integra- 
tion are given in Appendix \M Here, I remark that the code uses a space-centering 
discretization with the extensive quantities evaluated at the interfaces and the in- 
tensive quantities in the midpoints of the grid zones. Two time steps are adopted in 
each cycle: one to advance the velocity, and the other to advance the material state 
variables. The time step is chosen to comply with stability and accuracy constraints 
(see the discussion below). Special care is taken in the centering of the opacity 
when discretizing equation ( 13. 5 p in order to prevent numerical noise to appear due 
to the propa gation of the ra diation flow at the steep front where the opacity changes 
significantly (jChristyl . 119671 ). 



Depending on how the time derivatives of the differential equations are numer- 
ically evaluated, we have different schemes of integration. For example, given the 
equation. 



df{x,t) 

— Q-^ — = 9{f,x,t), (3.8) 



in a single time step the numerical scheme advances the solution from time t (where 
we known the solution) to time t + At. The advanced solution can in general be 
written as 
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fit + At) ^ fit) + At((l - 9) git) + 9git + At)). 



(3.9) 



If = the scheme is called explicit (or forward time differencing) and it involves 
only known values of the functions at time t. In other cases, with 6 between 0.5 and 
1, the scheme is called implicit and the solution involves the initial unknown function 
fit + At), requiring an iterative technique to find the solution. More specifically, 
when 6 = 1 the scheme is called fully implicit (or backward time differencing) and 
the case of 6' = 0.5 is usually known as centered time differencing. Although, the 
explicit schemes are easier to solve than the implicit ones, the former have a stringent 
requirement on the time step used in the calculation in order to ensure the stability 
of the numerical scheme. This requirement is known as the Courant condition., which 
states that At < Ax/v where Ax is the width of the grid zone and f is a characteristic 
speed which in hydrodynamics problems is adopted as the sound speed. Physically, 
this condition states that the information may not spread through the grid at a speed 
grater than the characteristic speed defined by the problem. 



In our code, we use an explicit scheme for the integration of the hydrodynamic 
eq uations, but a semi-im plicit scheme for the temperature, similar to the one used 
by 



Falk fc Arnett 



(119771 ). The equations are linearized in 6T and solved iteratively 
for each time step using the tridiagonal method (see Appendix |A]). The discretiza- 
tion typically uses 300 mesh points, with a finer sampling for the outer layers typ- 
ically smaller than 1O~^M0. This value was chosen based on tests performed with 



1988 



Ensman &: Burrows. 



19921 



our model and following previous works (iWoosley 
which show that the early light curve is sensitive to the mass zoning in the outer 
layers when a coarse grid is used. 



Given that we are using an explicit hydrodynamic scheme, the time step should 
be chosen as a fraction of the minimum of the Courant condition for all zones in 
order to achieve stability. I note that we have also imposed additional conditions on 
the time step: we have required that changes in temperature, density and flux over 
one time step be less than 5%. 
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The formation of the shock wave (SW) following core collapse is simulated 
by artificially adding internal energy ("thermal bomb") almost instantaneously in 
the inner region of the core. We usually set a mass cut of 1.4 Mq which the ma- 
terial assumed to collapse and form a neutron star or black hole. In order to do 
this, we have included an additional term in the energy equation 13.41 of the form, 
dEexp/dm = C exp{—t / t) exp{—m / mo) . This expression was arbitrarily chosen but 
its integral is exactly the explosion energy which is free parameter of our model. 
We have also tested explosions generated by injecting kinetic energy (in this case it 
is necessary to include an additional term in the Euler equation 13. 3p and we have 
obtained similar results. The latter method, however, leads to very short numerical 
time steps, and therefore to slower calculations. We have checked the accuracy of 
our calculation method by testing the conservation of energy. The total amount of 
energy is conserved within 0.6% but if we consider the conservation of energy be- 
tween two consecutive time steps, this is is within 4 xlO~^. We consider it very 
adequate for our purposes. 

3.1.2 Input physics 

The equation of state (EOS) is calculated using simple expressions for T, p, compo- 
sition, and for the ionization degrees of hydrogen and helium which are computed 
assuming local thermodynamical equilibrium (LTE). Therefore, the degree of ioniza- 
tion is determined by solving the corresponding set of Saha equations for ionization of 
hydrogen and the first and second ionization of helium. Ionization of heavy elements 
is neglected in our EOS (but, of course, it is taken into account in the calculation 
of the opacity). In general, our EOS consists of two terms, the radiation and gas 
contribution. The latter can be also divided into electron and ion terms. Thus, for 
example the total pressure is written as, 

P = PRad + Poas = PRad + Pe- + Pion (3.10) 

with 
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Rad 



4 R 
6 H 



(3.11) 



where fi is the molecular weight for the non-ionized material. The electron pressure, 
Pe-, is a much more complicated function. It is influenced by the degree of ionization 
(determined by Saha equation) and degenerac y effects. We calcu l ate th e pressure 
of degeneracy using the method described by iKippenhahn et al. although 



the condition of degeneracy is rarely reached in our calculation. We tested our 
results using a mor e sophisticated EOS, such as that of the Los Alamos Tables 
(IRogers et al.l . Il996l ) , and we did not find any significant difference with respect to 
the results obtained using our simple EOS. This is expected due to the low densities 
attained in most of the layers of the models. 



The Rosseland mean opacity, function of T, p (or P) and chemical 

composition. There are tables available in the literature that provide the opacity for 



a wide range of T, p and chemica 



comp osition. In our calculations we use the OPAL 



opacity tables (llglesias fc Rogersl . Il996l . and references therein). As these tables are 
giv en for T > 6 x 10^ K , we co mplement them with the opacity tables provided 
by [Alexander fc FergusonI (jl994[ ) for lower temperatures which includes molecular 
opacities. The tables are interpolated in order to guarantee a smooth transition at 
T = 10^ K from one to the other. 



These tables allow us to calculate opacities for several metallicities. Also, for a 
fixed metallicity, different mixtures of H, He, C and O can be used. Although in this 
work we adopt a fixed value of Z = 0.02, departures from this value (as expected in 
the inner regions of the object) are taken into account as excesses of C and O with 
respect to the values of the adopted metallicity, at the expense of He. 



Figure 13.11 shows the opacity given by the tables as a function of temperature 
for Z = 0.02 and different densities. The range of temperature shown corresponds to 
the values reached during the evolution of a SN II-P while the values of density are 
restricted to those achieved during the plateau phase. Also shown in the plot is the 
contribution to the opacity from electron scattering. Note that electron scattering 
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IS 



the dominant source of opacity for T > 10^ K and p < 10 g: 



—in — s 



The Rosseland mean opacity includes scattering and absorption processes. Scat- 
tering is the dominant process (see Figure I3.ip in the supernova ejecta until most 
of the electrons recombine with ions and absorption processes become an important 
source of opacity. On the other hand, in rapidly expanding envelopes where large 
velocity grad ients are present, the Rosseland mean opacity underestimates the true 
line opacity (IKarp et al.l . Il977l ). which hinders the estimation of the actual opacity 
in the outermost (recombined) layers. Another effect that is not included in the 
calculation of k is non-thermal excitation or ionization of atoms/ions which is pro- 
duced by Compton scattering of 7-rays emitted by radioactive decay of ^^Ni and 
^^Co. The LTE ionization used in the calculation of /t considerably underestimates 
the true ionization. The correct way to treat these effects is to calculate the actual 
contribution of non-thermal ionization to the opacity and to include the expansion 
opacity of lines. However, such treatment is beyond the scope of this study. We 
adopt an alternative approach that has be e n ext e nsively used in the l i teratur e (e.g. , 



Shigeyama fc Nomoto 



1990 



Herzig et al 



1990 



Swartz et al. 



1991 



Youngl . I2OO4J ) 



which consists of using a minimum value of the opacity (or "opacity floor"). The 
minimum opacity values adopted in this work are: 0.01 cm^ g~^ for the envelope 
material, and 0.24 cm^g~^ for the metal-rich core material. These values were cho- 
sen by performing a comparison with results of the STELLA code. We found an 
excellent overall agreement between the light curves g iven by our code and STELLA 



(IBlinnikov fc Bartunov 



1993 



Blinnikov et al. 



19981 ) both in terms of the duration 



of the plateau and the morphology of the bolometric LC) when the previous val- 
ues of the minimum opacity were adopted (see § 13.31 Note that STELLA is an 
implicit hydrodynamic code that incorporates multi-group radiative transfer, and 
additionally uses different o pacity tables and includes the effect of line opacities 
fisorokina fc Blinnikovl . I2OO2I ). 



26 




Figure 3.1: The run of the Rosseland opacity on temperature (T) and density (p) 
for solar metaUicity {Z = 0.02) as used in our calculation without including any 
"opacity floor" (see discussion in Section I3.1.2p . The ranges of T and p shown 
are typical for SNe II-P. The hydrogen (X) and helium (Y) mass fractions used are 
indicated. We also include the electron scattering opacity (considering full ionization, 
which is certainly unrealistic for the low temperature sector of this plot) in order to 
show the dominance of this source for T > 10^ K and p < 10~^° g cm~^. Note that 
this value of the density is reached early-on in the supernovae evolution. For lower 
temperatures the absorption processes become an important source of opacity. 
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3.1.3 Gamma-ray Deposition 



During the explosive nucleosynthesis produced in supernova explosions, unstable 
isotopes of iron elements are formed. The decay of these isotopes and subsequent 
thermalization of the decay products generate extra energy that contributes to power 
the LC. The most abundant radioactive isotope produced is ^^Ni which decays with a 
6.1 day half-life to ^^Co which in turn decays with a 77.7 day half-life to stable ^^Fe. 
The decays produce energetic 7-rays and positrons which are thermalized, providing 
the thermal luminosity of the SN. 



In order to include this source of energy in the calculations, it is necessary to 
determine a local heating rate which generally is not the same as the simple ra- 
dioactive decay rate. Instead, the decay rate should be modified by the probability 
of thermalization determined by the rate at which 7-rays and positrons deposit en- 
ergy at various points in the gas as they travel through the ejecta. This problem is 
very complicated if tackled from first principles. A set of 7-ray transfer equations 
must be solved simultaneously with the hydrodynamic equations which requires a 
multiple-energy group Monte Carlo calculation. However, if the ejecta is optically 
thick to 7-rays, it is possible to assume that the gamma rays deposit their energy 
locally. Several studies of SNe II-P have previously assumed this hypothesis consid- 
ering that the energy contribution of radioactive material is only relevant toward the 
end of the plateau phase. Such approach would be correct only if ^^Ni was deeply 
concentrated in the ejecta in which case gamma photons could hardly diffuse out 
of the regions where they form. We believe there is no just ification to assume this 



type of ^^Ni di s tribu t ion, as shown i n studies of SN 1987A (iShigeyama et al 



1988; 



Woosley et al. 



1988 



Arnett 



11988: 



Blinnikov et al 



2000 



among others) and the 
diffusion of gamma rays must be properly calculate from the location where they 
are emitted to the outer regions. To calculate this, we solve the gamma-ray trans- 
fer in the gray approximation for any spherically symmetric distribution of ^^Ni, 
assuming that gamma rays interact with matter only through absorption. Details 
of this calculation are presented in Appendix [Bl It has been shown by comparison 
with Monte Carlo simulations of Type la supernovae that the complex scattering 
process between gamma rays and electrons can be satisfactorily approximated as an 
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absorptive process (ISutherland &: Wheelerl . Il984j : ISwartz et al.l . Il995l ). The adopted 
value for the gamma-ray opacity is = 0.06 ?/e cm^g~^, where Ue is the number of 
electrons per baryon. 



The energy rate per gram released by Ni-Co-Fe decay is 



e,ad = 3.9 X 10^° exp(-t/rNi) + 6.78 x 10^[exp(-t/rco) - exp(-t/rNi)] erg g~S-\ 

(3.12) 

where = 8.8 days, and tqo = 113.6 days are the mean lifetimes of the radioactive 
isotopes. The amount of energy deposited at each point is given by the solution of 
the gamma-ray transfer multiplied by the previous expression. In Figure [32] we show 
the gamma-ray depositior0 profile for the case of a polytrope with index n = 3, initial 
mass 10 Mq and different initial radii. The left panel is for a constant distribution 
of ^^Ni up to 3Mq and the right panel is for an exponential distribution of ^^Ni, also 
up to 3Mq. Note that the diffusion of gamma rays from the region where they form 
becomes more noticeable as the object becomes more extended and diluted. 

The possibility of using an arbitrary distribution of ^^Ni allows us to study 
different types of mixing and their effect on the resulting LC. As shown in 15.1.31 
we find that radioactivity becomes an important source of energy, even during the 
plateau phase, if we allow an extensive mixing of ^^Ni. 



3.1.4 Initial models 



There are two different ways to determine the initial (or pre-SN) models: those 
coming from stellar evolution calculations ("evolutionary" models), or those from 
non-evolutionary calculations where the initial density and chemical composition 

^The deposition of gamma rays is defined by the energy deposited at each point normahzed to 
the value corresponding to complete thermalization at the same location where the gamma rays 
are emitted. 
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Figure 3.2: Gamma-ray deposition as a function of mass for a polytrope with index 
n = 3, initial mass 10 Mq, and different initial radii. (Left) for a constant distribu- 
tion of ^^Ni up to SMq] (right) for an exponential distribution of ^^Ni, also up to 
3Mq. The diffusion of the gamma rays out of the region were they form becomes 
more noticeable as the object gets more extended and diluted. 

are parameterized in a convenient way. Both types of models have been tested with 
our code. 

As non-evolutionary models, we have calculated single and double polytropes. 
The single polytropic models were numerically computed by solving the "Lane- 
Emden" equation for any polytropic index (n). This allowed us to study many 
different configurations for different values of n, initial masses and radii, in hydro- 
static equilibrium. Some examples of initial density profiles obtained using different 
polytropic indices are shown in Figure [3731 

Although a single polytrope may represent very well the envelope of the real 
pre-supernova model, the inner dense part which is expected for this type of objects 
is not well reproduced with this simple initial model. One way to improve this 
situation is to consider two polytropes: one representing the inner dense core, and 
the other accounting for the outer extended envelope. 
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Figure 3.3: Density distribution as a function of mass for a single polytropic model 
using different polytropic indices, ra, which are indicated in the plot. The total mass 
and radius in all the models are 19 Mq and 800 Rq. Note that higher values of n 
produce more compact configurations. 



The numerical calculation of double polytropic models is not trivial because the 
solution is extremely sensitive on the adopted initial guess value of the free param- 
eters. Specifically, it is necessary to solve the equations of hydrostatic equilibrium 
and mass conservation assuming a polytropic approximation where the pressure has 
a dependence on density of the form: P = Kp^ = Kp"/("+i), where is a constant 
and n is the polytropic index. Two different polytropic relations with different in- 
dices (n, and Uo) and constants {Ki and Ko) are adopted to mimic the characteristic 
structure of a red supergiant. To calculate the composite polytrope, the previous 
equations are numerically integrated outward from the inner border (stellar core) to 
a pre-selected point (fitting point) using the internal polytropic relation. A second 
inward integration is done from the outer border (stellar surface) to the fitting point 
using the external polytropic relation. The problem is equivalent to a two-point 
boundary condition for the case where there are unknown free parameters at both 
ends of the domain. In our case, the free parameters are rij, and Ko, while the 
external index is fixed to no = 3. The inner boundary conditions are the central 
density (pc) and the core mass (m^), while the external boundary conditions are the 
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initial mass (M) and the external density (p^). The integrations are performed along 
the radial variable r from to R (radius of the object). The values of r^, -R, M and 
Pc are fixed for each configuration, while rric, Po are functions of the free parameters 
of the problem. Starting from initial guesses, the values of the free parameters are it- 
eratively sought so that the solution joins smoothly at the fitting point. This process 
is performed using a shooting method (see Numerical Recipes, chapter 17). Once 
the values of rij, Ki and Ko are found, the density (or pressure) distribution and 
mass distribution can be calculated. The initial temperature profile is calculated in 
an iterative fashion using our EOS to ensure hydrostatic equilibrium. This prevents 
the formation of a spurious shock wave. 

Using this method a variety of initial models in hydrostatic equilibrium can be 
calculated by changing the fitting point, the initial mass, the radius or the central 
density as shown in Figures 13.41 and 13.51 The effect of the initial mass and fitting 
point on the density distribution is shown in Figure 13.41 The effect of the initial 
radius is shown in Figure 13.51 Note that different configurations can be obtained by 
changing the fitting point even if the initial mass and radius remain fixed. The same 
is valid if one changes the central density or the external polytropic index. 

These parametric profiles allow us to study how the internal structure of the 
progenitor affects the LC. Figure 13.61 shows a comparison between LCs obtained 
using a single and a double polytropic model for an initial mass of 19 Mq and an 
initial radius of 800 Rq. During the early part of the LC evolution {t < 20 days) 
both models are very similar. Significant differences appear later on. Note that the 
single polytropic models produce a prominent bump before the sharp decrease in 
luminosity which marks the end of the plateau. This bump is a consequence of the 
lack of a dense core in the initial model and is never observed in SNe II-P. 



We have also calculated LCs using initial models derived from stellar evolution 
calculations. Pre-supernove models of four different sources were used: (1) "Su - 



perno va Science Center-UCSC" (http://www.supersci.org/ucsc/), (2) iLimongi et al. 



(|2000[ ) (http://web.oa-roma.inaf.it/localinfo/staff/webhost/limor igi/data.html), (3) 
[Com ment: Pedir cita a Omar:Benvenutto et al.] and (4) 



Umeda fc Nomoto 



(120051). Figure 13171 shows the density distribution for some of these models. Also 
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m [Mq] m [Mq] 

Figure 3.4: Density distribution of double polytropic models as a function of mass 
for a model with radius of 800 Rq. (Left) The effect on the initial density of varying 
the mass. The values used in the calculation are indicated in the plot. (Right) The 
effect on the initial density of varying the fitting point for an mass of 19 Mq. Note 
that even if mass and radius are the same, different configurations are obtained by 
changing the fitting point. 
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Figure 3.5: Density distribution of double polytropic models as a function of radius. 
The different curves represent different values of the radius used in the calculation 
for a fixed mass of 19 Mq, as indicated in the plot. 
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Figure 3.6: Comparison between the LC obtain using a single (red line) and double 
(blue line) polytrope. While the early LCs {t < 20 days) are very similar, important 
differences appear toward the end of the plateau. 



shown for comparison the density distribution of a double polytropic model. The 
LC obtained using a pre-SN model of source (1) is shown in Figure [331 Note that 
some bumps appear around t = 35 days in the LC which are not present in observa- 
tions. In § 13.31 a discussion using a pre-SN of source (4) is presented. 



Having tested the code with various evolutionary initial models, I generally find 
better agreement with observations when using our parameterized double politropic 
profiles, in concordance with previous studies that emp l oy similar methods to com- 



pute the initial structure (jUtrobin 



1993 



Baklanov et al. 



2005 



Utrobi 



2 



20071 . among 



others). Thus is the approach adopted throughout the rest of this thesis. 
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Figure 3.7: Initial density distributions form stellar evolution. The pre-SN models 
come from different sources, as indicated the text. A double polytropic model is also 
included for comparison. 
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Figure 3.8: LC calculated using the initial model from stellar evolution calculations 
of the "Supernova Science Center-UCSC" . Note that some bumps appear in the LC 
around t = 35 days which are not present in observations. 
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3.2 Code structure 



The code is organized in several subroutines with different hierarchy. Figure 13.91 
schematically shows how our computational program is organized. 



SETUP 

M, R, E, Mm, 

to, tend, AtQ,... 



MAIN DRIVER 

t= t + At 



if t < tend 



at some t 



OUTPUT 

r{m), p{m), T{m), v{in) 



iit = tn 



INITIAL 
MODEL 



Single Polytrope 
Double Polytrope 
Evolutionary 



PHYSICS 
DRIVER 



^-RAY TRANSFER 



TIME STEP 

At 



HYDRODYNAMICS 



TRANSPORT 



EOS 



OPAL 



Figure 3.9: Code scheme 



The main routine or "MAIN DRIVER" controls the program during the course 
of a calculation through a series of subroutine calls which are governed by logical 
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sentences based in part on the input information. The input of the program is rep- 
resented in the "SETUP" block and includes the information necessary to generate 
a specific calculation. There, an input file is read containing the following informa- 
tion: type of initial model to use in the calculation (polytrope, double polytrope or 
evolutionary model), initial mass (M), radius (R), nickel mass (Mni), distribution 
of ^^Ni and explosion energy (E). The initial time step (Ato) and the ending time 
for the model calculations {tend) are also included in this file. Other parameters 
which control the accuracy of the solution and the number of iterations allowed for 
the convergence of the iterative transfer calculations are also included there. The 
"PHYSICS DRIVER" block represents the piece of the code which calls subroutines 
that calculate the evolution of the physical quantities. Specifically, the following 
subroutines are called from this driver: (1) a routine to calculate the explicit hydro- 
dynamic equations, (2) a transport routine which calculates the diffusion radiative 
transfer using an iterative method as described in Appendix [Xj and (3) the routine 
for solving the gamma-ray transfer (see Appendix [B]). The EOS and OPAL sub- 
routines are called from the transport routine and represent the calculation of the 
equation of state and opacity. The results of the calculations, i.e., the model struc- 
ture, r(m), P(m), p{m), T{m), etc., and the evolution of certain quantities such as 
bolometric luminosity, Lhoi(t), effective temperature, T^fiit), photospheric velocity, 
Vph{t), among others are written in output files at specific moments during the evo- 
lution. The writing process is schematically represented by the block "OUTPUT". 
This subroutine is called from the "Main Driver" at the specific moments defined 
in the "SETUP" routine. After each step of calculation, a new step of time. At, is 
chosen in the "TIME STEP" block depending on stability and accuracy conditions. 
The "MAIN DRIVER" routine determines whether a new calculation is necessary, 
basically by checking that the new time step, t = t + At is less than tend- In that 
case, the whole process is repeated for time t + At. 

Figure 13.101 shows a typical bolometric LC obtained with this code. The initial 
model used in the calculation is a double polytrope with initial mass 19 Mq, initial 
radius 800 Rq and injected energy 1 foe (1 xlO^^ erg s~^). From the Figure it 
is possible to distinguish various phases: (1) a maximum in luminosity or " shock 
breakout" , (2) a phase where the luminosity is nearly constant or "plateau" , (3) a 
transition phase, and (4) a radioactive tail. These phases are studied in detail in 
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Chapter O 




Figure 3.10: A typical bolometric LC of a SN II-P. The initial model used in the 
calculation was a double polytrope with initial mass 19 M©, initial radius 800 Rq, 
and injected energy 1 foe. The different phases during the LC evolution are indicated 
in the plot. 



3.3 Discussion on the approximations 



Several approximations are made in the equations of radiation hydrodynamics. First, 
we assume that the fluid motion can be described by one-dimensional, radially sym- 
metric flow. The explosion mechanism of core-collapse SNe is not well known. It 
may be an very asymmetric process. However, for this particular subtype of super- 
novae with very extended hydrogen envelopes the asymmetries expected from the 
explosion mechanisms itself appears to be smoothed. T his is supported by recent 
spectropolarimetric studies ( iLeonard fc Filippenkd . l2005l ). 



We use the equilibrium diffusion approximation to describe the radiative trans- 
fer. This approximation assumes that radiation and matter are strongly coupled 
with a single characteristic temperature and a spectral energy distribution described 
by a black body function (BB hereafter). The approximation breaks down at shock 
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breakout and at late phases when the ejecta is completely recombined and the ob- 
ject becomes transparent. Fortunately, during the plateau phase — which is our 
main interest — , this approximation is fairly adequate. Also note that non-LTE cal- 
culations of SN II spectra show that deviations from LTE have significant effects on 



spect ral lines but not on the overall continuum ( iBaron et al 



1996; 



Dessart fc Hillierl . 



20081 ). At late phases, for SNe that experience little interaction with the interstellar 



medium, as is the case of SNe II-P, the bolometric luminosity can be simply ap- 
proximated as the luminosity deposited by the radioactive decay of ^^Co, at least 
during the early part of the radioactive tail. This is supported by the fact that the 
luminosity declines obeying an exponential law with a very similar rate to that of the 
decay of ^^Co. The radioactive scenario has bee n directly confirn ied for SN 1987A 
through the detection of 7-ray lines from ^^Co (IMatz et al.l . Il988[ ). No attempt to 
apply our model beyond ~180 days is done. 



During the transition between plateau and radioactive tail, the envelope is fully 
recombined and the notion of a photosphere loses meaning. Thus, this transition 
regime is poorly described with our radioactive transfer prescription and therefore 
a detailed analysis of this phase cannot be asessed here. As stated above, we have 
performed a comparison of our calc ulations with those of th e STELLA code using 
the same initial model provided by lUmeda &: Nomotd ( 120051 ). We obtained an ex- 
cellent agreement between the bolometric luminosities derived by both methods, as 
shown in Figure 13.111 This is very satisfactory considering that the STELLA code 
involves a more sophisticated treatment of the radiative transfer by solving these 
equations using a multi - group prescription and including t he effect of line opacities 
( Blinnikov fc Bartunovl . Il993l : ISorokina fc Blinnikovl . I2OO2I ) . 



Finally, I remark that the very good agreement obtained between both codes 
was reached using the values of opacity minimum (/tmin) mentioned in § 13.1.21 i.e. 
for the envelope material Kmin = 0.01 cm^ g~^. For other values of the opacity min- 
imum commonly used in literature, the comparison with the STELLA code, which 
incorporates the effect of expansion opacity in the correct way, was not satisfactory. 
In order to show the effect on the LC of the value adopted for the opacity floor, I 
calculated models using three different values for the envelope material: Kmin[cm^ 
g~^] = 0.001, 0.01 and 0.1, as shown in Figure 13.121 Clearly, from this Figure, there 
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Figure 3.11: Comparison between the LC obta ined with our code and th e STELLA 
code. The initial model used was provided by lUmeda fc Nomotd ( l2005l ). Note the 
very good agreement between both LCs despise the simplifications employed here. 
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is a critical effect of the opacity minimum on the resulting LC. Hence, it should not 
be taken as a free parameter. 
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Figure 3.12: LC for three different values of the opacity minimum {Kjnm, expressed 
in cm^ g~^ in the plot). Note the critical effect of the adopted opacity minimum 
on the LC. The value finally used in our code, Kmin = 0.01 cm^ g~\ was chosen by 
comparison to the STELLA code (see Figure 13. lip 
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Chapter 4 

Comparison with observables 



Since our code produces bolometric light curves and effective temperatures (Te//) it 
proves necessary to calculate these quantities from the observed photometry in order 
to compare our model with observations. Although there are many SNe with optical 
observations, only a handful objects are available with observations over UV, optical 
and IR wavelengths. The purpose of this section is to use those SNe with observed 
bolometric light curves and explore the feasibility to derive a bolometric correction 
(BC) for other SNe with optical observations alone. 

In §4.H I present the calibrations derived for BC and Te// from BVI photometry 
using data of three well-observed SNe II-P, SNe 1987A, 1999em and 2003hn, and two 
sets of atmosphere models by Eastman et al. (1996) and Dessart & Hillier (2005b) 
(E96 and DOS hereafter, respectively). The typical scatter of the BC is 0.11 mag 
and allows us to calculate bolometric luminosities for many other SNe II-P having 
BVI photometry alone, opening thus the possibil ity for a statistical analys is of the 
physical properties of this type of object (see also lBersten &: Hamuyl ( 120091 )). 
Additionally, the photospheric velocity is another useful parameter to compare with 



observations. In § 14.21 I describe the prescription used to compare the observed 
velocities with the photospheric velocities provided by our models. 
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4.1 Analysis of BC and Tg// -color relations 



I begin this section by describing the observational and theoretical material used 
to calculate BC and Teff- Then, in § 14.1.11 1 proceed to calculate the bolometric 
luminosities for these data. Finally, in § 14.1.21 and § 14.1.31 1 derive calibrations for 
BC and Teff as a function of colors. Along this section, I refer to the SN evolution 
in terms of time or color indistinctly. This is well justified during the plateau phase 
in which the SN atmosphere expands, cools and monotonically turns redder. 



In order to examine if a bolometric correction can be derived from optical colors, 
we made use of the three SNe II that possess the best wavelength and temporal 
coverage. Two of these are genuine SNe II-P, SN 1999em and SN 2003hn, and the 
third is the famous SN 1987A which, except for its peculiar light curve, shares most 
of the spectroscopic properties of SNe II-P. Most of the data of these three SNe 
were obtained at Cerro Tololo Inter- American Observatory (CTIO), Las Campanas 



Observatory (L CO), and the European Sou ther n Observatory 



more details see 



Hamuy fc Suntzefla ( 1l990l ) and 



Hamuv et al.l (j20l3) for SN 1999em, and 



Bouchet et al 



Krisciunas et al 



ESQ) at La Silla. For 
fll989ll for SN 1987A, 
for SN 2003hn. 



The photometric bands used in this analysis were UBVRIZJHK for SN 1999em 
and SN 2003hn, and UBVRIZJHKLM for SN 1987A. 



We adopt C epheid distances f or SN 1987A and SN 1999em with correspo nding 
values of 50 kpc ( Freedman et al. . 2001 ) and 11.7 Mpc fjheonard et al. . 2003). For 
SN 2003hn we used a distance of 16.8 Mpc, as derived by Olivares et al. 2008 using 
the Standardized Candle Method. We correct the photometry for Galactic and host- 
galaxy extinction. We perform such corrections using Galactic visual absorptions 
of Ag^^ =0.249 for SN 1 987A , A^^-^=0.13 for SN 1999em, and A^^^=0.043 for SN 



2003hn ( ISchlegel et al.l . 119981 ) . assuming a standard reddening law with Ry= 3.1 
as given by Cardelli et al. (1989). The values used f or host-galaxy absorption are. 



A host_ 

I\y - 



0.216 for SN 1987A, A^r *=0.18 for SN 1999em ( ]Hamuvl . l200lh and A 



for SN 2003hn ( jPessartl . l2008[ l . again assuming Rv= 3.1. 



,^«^t=0.56 
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We also used in this analysis spectral energy distributions (SEDs) from two sets 
of SN atmosphere models (E96 and DOS). These models depend on several param- 
eters such as, luminosity, density structure, velocity and compositio n. For more de- 
tails o n the input parameters of such models the reader is referred to lEastman et al. 



and bessart fc Hilhej ko05ii \. We use a total of 61 model spectra from E96, 
and 107 from DOS. We discard 31 spectra from DOS which do not have enough UV 
coverage. 



4.1.1 Bolometric Luminosity Calculations 



By definition, the bolometric luminosity is the integral of the flux over all frequencies. 
This integration can be done in a straight-forward way for the spectral models of E96 
and DOS summing the flux over wavelength. With the purpose to estimate bolometric 
corrections and colors for the models we compute BVI s ynthetic magnitudes using 
the filter transmission functions and zero points given by iHamuyl (l200ll ). 



For the three well-observed SNe the calculation of bolometric luminosities is 
performed from reddening-corrected broadband magnitudes using the values men- 
tioned in section W7L[ K-corrections are neglected due to the small redshifts involved. 
We began by computing a quasi-bolometric light curve using all the available broad- 
band data. The magnitudes were converted to monochromatic fluxes at the specific 
effective wavelength of e ach filter using the transmission functions and zero points of 
the photometric system (iHamuyl . l200ll ). At epochs when a certain filter observation 
was not available, we interpolated its magnitude in time using the closest points. 
The total "quasi-bolometric" flux, Fg^oi was computed using a trapezium integration 
over U -K for SN 1999em and SN 2003hn, and over [/ - M for SN 1987A. 



To estimate the missing flux in the UV and IR, Fuy and Fjji, we fit at each 



epoch a blackbody (BB) function to the monochromatic fiuxe4j as shown in Fig- 
ure 14.11 for the case of SN 1999em at ~ 7 days post explosion. At early epochs the 



^ These fits are restricted to the plateau phase where the envelope of the SN was optically thick, 
and also to the transition to the nebular phase. On the radioactive tail we did not use any UV or 
IR corrections. 
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Figure 4.1: Blackbody (BB) fit to the monochromatic fluxes for the case of 
SN 1999em at ~ 7 days after explosion. The BB function provides a very good 
representation of the flux. The Tbb and 6 shown in the inset are the color temper- 
ature and angular radius of the SN 1999em yielded by the BB flt. 

BB model provided very good flts to the fluxes in all bands. As the photosphere 
became cooler the f/-band flux started to depart from the BB model in which cases 
we exclude this point from the flt. At later epochs, subsequently the S-band and 
l^-band data points show the same behavior, departing from the BB model. The 
reason for this is related to the strong line blanketing that develops with time in that 
part of the spectrum. 

On the IR side the flux is extrapolated to A = oo using the BB flts described 
above. The integral of that function between the longest observed effective wave- 
length and A = oo is adopted as the IR correction (Fjr). This correction increase 
with time but always remains below 7% for the three SNe. 

On the UV side, we extrapolate from the effective wavelength of the U band 
to A=0 using the BB flt on all epochs except when the ?7-band flux falls below the 
BB model. In these cases, we extrapolate the [/-band flux using a straight line to 
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Figure 4.2: Contributions to the bolometric flux as employed this work. Fu^x is 
the "quasi-bolometric" flux computed by integration of all broad band available by 
each SN. Fuv and Fjji represents the UV and IR contribution calculated using a 
extrapolation of the BB fit to the broad band or a linear extrapolation. 

zero flux at 2000 A. Our choice of A= 2000 A as the wavelength where the flux 
goes to zero was based on the behavior of the atmospheric models for which the flux 
blueward of 2000 A is negligible in comparison with the total flux. The integrated 
flux under the Planck function (or straight line) between the effective wavelength of 
the U fllter and A=0 (or A=2000 A) is taken as the UV correction (Fuv)- Figure 
shows schematically the different contributions to the bolometric flux employed in 
this work. 

The size of the Fuv correction relative to the total flux for the three SNe and 
the two sets of atmospheric models is shown in Figure 14.31 as a function of {B — V). 
The flrst thing to note is the overall good agreement in F^y between the observed 
SNe and the atmospheric models. Second, it is evident that the UV correction is 
very important at early epochs and becomes nearly irrelevant at the latest epochs. 
Thirdly, note that in the very blue end, where the UV corrections are of order 50- 
80%, there is some disagreement between the atmosphere models and SN 1999em. 
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Figure 4.3: UV contribution to the total flux as a function of {B — V) for SN 1987A, 
SN 1999em and SN 2003hn and for the models of E96 and DOS. At early times, 
when [B — V) ^ 0.2, the UV flux represents a significant fraction of the total flux, 
implying larger uncertainties in the UV correction. 

The spectral models suggest larger UV corrections and, as argued below, they are 
more trustworthy at these early epochs than the extrapolation of the broadband 
magnitudes. To prove this point, we calculate Fjjv for the models using the same 
technique that we employe for SN 1999em, i.e. by computing synthetic magnitudes 
for all passbands, converting them to monochromatic fluxes, and fitting a BB to the 
resulting points (instead of the direct integration of the SED. In this case the UV 
correction proves closer to the UV correction derived from SN 1999em. We conclude 
that the UV extrapolation using the BB fits to the broadband magnitudes at very 
early epochs underestimates somewhat Fuv, so we end up using only the atmosphere 
models at such epochs. At later times {B — V> —0.04), where the differences between 
data and theory become negligible, we adopt both the models and the observed data. 

The sum Fgboi+Fuv+Fm yields the bolometic flux Fboi (the integration under 
the solid curve of Figure 14. 2p . Then we transform flux into luminosity using the 
distances given in section 14.11 The resulting bolometric luminosities for SN 1987A, 
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SN 1999em and SN 2003hn are shown in Figure I4.4[ As a comparison, the sohd hne 
shows the bolometric luminosity of the SN 1987A obtained by Suntzeff & Bouchet 
(1990). I find very good quahtative agreement between both bolometric light curves 
for SN 1987A. There is a systematic difference which remains smaller than 0.04 
dex at all times between both calculations. Such differences are consistent with the 
uncertainties arising from the use of different photometric data sets and different 
integration and interpolation scheme. 



As can be seen in Figure 14.41 the morphologies of the bolometric light curves 
for the three SNe are very different, especially that of SN 1987A which shows a 
broad maximum, not observed in classical SNe II-P, and a less luminous light curve 
(up to the transition to the radioactive tail). The peculiar light curve of SN 1987A 
is well known and is attributed to the fact that its progenitor was a compact blue 
supergiant that l ead to a dim in i tial plateau and to a light cu rve promptly powered 
by radioactivity (IWoosleyl . Il988t IShigeyama fc Nomotd . [1990) ■ However, the three 



objects showed a similar initial phase of rapid fading and cooling until the outermost 
parts of the ejecta reached the temperature of hydrogen recombination (adiabatic 
cooling phase). A second phase can be distinguished for SN 1999em and SN 2003hn 
which corresponds to the plateau where the luminosity remained nearly constant 
while hydrogen was recombining. The duration and the slope of the light curve 
during this phase were different for each supernova. This is related to the properties 
of the progenitor object, mainly to the mass and radius of the hydrogen envelope. 
The shape of the light curve for SN 1987A during this phase was very different as 
mentioned above. It showed a broad maximum characterized by a slow rise of ~ 90 
days followed by a more rapid decline for about 30 days. Finally, I can distinguish a 
third phase, the radioactive tail which has similar slope for all three SNe. Here, the 
luminosity exhibits a linear decline, and it is dominated by the radioactive decay of 
^^Co into ^^Fe . The luminosity in this part of the light curve is a direct indicator of 
the amount of ^^Ni synthesized in the explosion (jWoosley et al.l . Il989l ). the parent 
product of ^^Co. We deduce from this that SN 1987A produced more ^^Ni than the 
other two SNe. 



48 




Figure 4.4: Bolometric luminosity of SN 1987A, SN 1999em and SN 2003hn com- 
puted from the integration of broadband optical and near-infrared data plus UV 
and IR contributions as explained in section I4.1.1I For comparison we include the 
bolometric luminosity of SN 1987A obtained by Suntzeff & Bouchet (1990) (solid 
line) . 
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4.1.2 Bolometric corrections versus Color 



There are many SNe that lack IR and UV observations for which it is not possible 
to calculate the bolometric luminosity using the method described in the previous 
section. For these cases it is necessary to know the bolometric correction required 
to convert a l^-band magnitude into a bolometric flux, i.e., 



BC = m,oi - [V - Avl (4.1) 

where Ay is the total visual extinction and mi,oi is the bolometric magnitude in 
the Vega system. Note that, since BC is deflned as a magnitude difference, it is 
independent of the distance assumed for each object. 

We calculate BC at all epochs for each of the calibrating SNe and all of the SN 
models using the bolometric luminosities computed in section 14.1.11 The bolometric 
fluxes were converted into Vega magnitudes in the following manner. 



niboi = -2.5 logio F^oi + 11.64, (4.2) 

where the zero point is obtained by integrating the SED of Vega given by Hamuy(2001) 
and forcing the resulting magnitude to vanish, i.e., mf,oz(Vega)= 0. 

We analyzed the dependence of the BC on color, using {B — V), {V — I) and 
{B — I). The main reason why w e did not try colors involving the R band is because 



our sample (iHamuy et al.l . |2010| ) has few SNe with /2-band observations. Figures 
14.51 and 14.61 show the resulting corrections as a function of the mentioned colors 
(corrected for dust) for SN 1987A, SN 1999em, SN 2003hn, and the models of E96 
and D05. In each of these plots the vertical bar indicates the approximate color 
corresponding to the end of the plateau. To the red of this mark are shown the BC 
corresponding to the transition between the plateau and the radioactive tail for the 
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Figure 4.5: Bolometric corrections versus {B — V) (left) and {V ~ I) (right) for SN 
1999em (open circles), SN 2003hn (stars), SN 1987A (filled circles) and the models 
of E96 (triangles) and D05 (squares). The vertical lines indicate the color at the end 
of the plateau phase. The solids line show polynomial fits to the points. The dashed 
curves correspond to the bolometric corrections of a blackbody spectrum. 



three SNe (no atmosphere models cover this phase). We do not include any of the 
nebular data in these diagrams, since the BB fits are not appropriate to extrapolate 
UV or IR fluxes at these epochs. 

These plots reveal a remarkable correlation between BC and intrinsic color, 
both for the objects and the models. It is very satisfactory that, even though SN 
1987A has a very different light curve compared with normal SNe II-P, it matches 
quite well the behavior of the other two SNe and the models. At very early times the 
BCs are quite large owing to the relatively larger flux contribution of the UV. During 
most of the plateau the BC remains very small around a value of zero, which implies 
that the V magnitude provides a very close proxy bolometric magnitude. During the 
transition from the plateau to the radioactive tail (redward from the vertical bar) 
the BC starts to depart from zero due to the larger flux contribution in the IR. At 
this phase SN 1987A shows some discrepancies, at the level of ~0. 1-0.2 mag, with 
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Figure 4.6: Bolometric corrections versus [B — I) for SN 1999em (open circles), SN 
2003hn (stars) , SN 1987A (filled circles) and the models of E96 (triangles) and DOS 
(squares). The vertical line indicates the color at the end of the the plateau phase. 
The solid line shows a polynomial fit to the points. The dashed curve corresponds 
to the bolometric corrections of a blackbody spectrum. 
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respect to the other two SNe. Since the bolometric fluxes for SN 1987A comprise 
two more IR bands, L and M, we examine the possibility that these discrepancies 
could be due to this fact. For this, we exclude the L and M bands for SN 1987A and 
recompute the bolometric flux in the same manner as for the other two SNe, i.e, by 
calculating Fm as the extrapolation of a BB fit to f/ — i^' photometry. This exercise 
shows that, while the BC corrections over the plateau phase do not change in any 
significant way (lending support to the Fm derived from BB fits), by the end of the 
plateau and at later times the new BCs increase and get closer to the other two SNe. 
The conclusion is that the differences observed during the transition are due to an 
inaccurate estimate of Fi^ from the BB fits restricted to U — K photometry. Adding 
L and M photometry at these late epochs does help and provides a more accurate 
estimate of Fjr. Therefore, during the transition we decided to exclude the BCs 
derived from SN 1999em and 2003hn. 

A good representation of the correlation between BC and colors can be obtained 
with polynomial fits of the form, 

n 

BC{color) = ^ai {color)\ (4.3) 

i=0 

where the order n varies for each color. Table 14.11 lists the coefficients obtained for 
the fit of each color, their range of validity and the number of data points used. The 
fits have dispersions (rms) of 0.11 mag for [B — V), 0.11 mag for {V — I) and 0.09 
mag for {B — I) in the whole range (plateau plus transition to the radioactive tail). 

The corresponding polynomials fits are also shown with solid lines in Figure 14.51 
and 14.61 As a comparison, the dashed lines in these Figures show the BC derived for 
a blackbody. The blackbody models represent well the data at early times (bluest 
colors), but evidently differ from the atmosphere models and the observed SNe at 
later epochs. 

As argued in section 14.1.11 we have good reasons to trust more the atmosphere 
models than the early data of SN 1999em, so we decided to exclude the latter from 



53 



Table 4.1: Coefficients of the fits to BC(color)^ 



B-V V -I B-I 



ao 

as 
a,4 
as 
ae 
ay 
as 



-0.823 -1.355 -1.096 

5.027 6.262 3.038 

-13.409 -2.676 -2.246 

20.133 -22.973 -0.497 

-18.096 35.542 0.7078 

9.084 -15.340 0.576 

-1.950 ••• -0.713 



0.239 
-0.027 



ranges 
No. points 
rms [mag] 




0.2,1.65] [-0.1,1] [-0.4,3] 
512 465 512 

0.113 0.109 0.091 



a 



BC{color) = J2 



n 



{colorY 



i=0 



our fits. Therefore, our cahbration should be considered more uncertain here. A 
further complication at early phases is the steep dependence of the BC on color. 
This means that a slight error in the measurement of the color, such as that due to 
a poor extinction determination, could cause a significant error in the determination 
of the BC. This problem is less pronounced if we use {V — I) to estimate BC at these 
epochs. We also test if a bolometric correction with respect to the R band instead 
of V would improve the situation, but we do not find any improvements. 

Using the coefficients given in Table 14. H it is possible to derive a bolometric 
luminosity for any SN II-P using only two (or three in the case of the {B — I) 
color) optical filters. If one knows the extinction and the distance to the object the 
bolometric luminosity can be computed as follows: 



logL[erg s-i] = -0.4 [BC (color) + V- Aotai{V) - 11.64] + \og{47iD^), (4.4) 



where D is the distance in cm to the SN and AtotaiiV) is the total, host plus Galactic 



54 



visual extinction. Note that combining this equation with equations ( 14. ip and ( I4.2p . 
the luminosity becomes independent of the arbitrary zero points chosen for the Vega 
magnitude scale, and it only depends on the observed flux density, color, extinction 
and distance. 

The calibrations of BC versus colors shown above are only valid during the 
optically thick phases since they involve BB fits to the photometry. In the nebular 
phase, we calculated BCs for the three SNe using the integrated flux between the 
observed bands (i.e. U - K for SN 1999em and SN 2003hn, and U - M for SN 
1987A). We did not attempt to add any flux beyond these limits since we did not 
have any physical model to extrapolate. As shown in the left panel of Figure 1^771 the 
BC for SN 1987A is almost independent of color, with a value ~ —0.7 mag and a 
scatter of only 0.015 mag. The other two SNe yield BCs 0.2-0.3 mag higher, with 
a slight dependence on color. We investigate whether these differences could be due 
to the inclusion of the two additional bands for SN 1987A: we removed the L and M 
bands from the BC and, not surprisingly, the agreement proved much better (Right 
panel of Figure I^Tll . We conclude that the L and M contributions to the bolometric 
flux is not negligible at the nebular phase. Hence, we take the value of —0.70 derived 
from SN 1987A as the best estimate of the BC at the onset of the nebular phase. 

We conclude this section with the claim that we have implemented a robust 
method to estimate BCs for SNe II-P which allows one to derive bolometric lu- 
minosities. If we trust the late behavior of SN 1987A as being representative of 
SNe II-P in general, our analysis implies an overall accuracy of 0.05 dex in BCs. 
Clearly, it would be interesting to check this result using L and M photometric data 
of other SNe II-P, but such data are currently unavailable. Our calibrations have the 
potential to be applied to many SNe observed over a limited wavelength range. 

4.1.3 Effective Temperature- Color Relation 

Along with bolometric luminosity, the effective temperature is a critical parameter in 
the comparison of observations with hydrodynamical models. Each model spectrum 
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Figure 4.7: ( Left) Bolometric corrections obtained during the radioactive tail phase 
versus {B — V), as derived using U through M photometry for SN 1987A (filled 
circles), and U through K photometry for SN 1999em (open circles) and SN 2003hn 
(stars). (Right) Same as before, but excluding the L and M bands for SN 1987A. 
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of E96 and DOS has an associated effective temperature (Te//), defined by the relation 
L = AuRpfj^aT^j-j where L is the input luminosity of the atmospheric models and Rph, 
the photospheric radius, is an output of the models. I examined the dependence of 
Te// on {B — V) and (V — I) colors derived via synthetic photometry from the model 
spectra, as described in section I4.1.1I The purpose of this analysis was to provide 
a calibration between temperature and color which could be used to easily derive 
estimates of Te// from observed colors for any SN II-P. Note, however, that Teff does 
not have a direct physical meaning for this type of object. It is simply a convenient 
contact point between hydrodynamical models and observations. 



Figure US] shows the effective temperature versus synthetic {B — V) and {V — I) 
colors for E96 and DOS models. As expected, there is a tight correlation between 
these quantities for each set of models. At early epochs, when {B — V) ^ 0.2 and 
{V — I) < 0.3, both models show consistent values of the effective temperature 
within their internal dispersion. Later on, however, when the plateau phase is well 
established, there are systematic differences in the behavior of both sets, with the 
models of DOS giving larger effective temperatures. Similar differences have been 
reported in the literature with regard to the dilution facto rs calculated from both 
sets of models ( iDessart fc Hillierl . 



200Sb 



Jones et al. 



20091 ). but there has been no 
clear explanation for the discrepancies. Note that Tg// ( oc ^jRph) is an output of the 
atmosphere model and depends on complicated details of the solution of radiation 
transport through the envelope, such as non-LTE treatment of the different species 
and metal line opacities. 



In order to represent the correlation between Tg// and colors shown in Figure 
14. 8[ we fit polynomial functions of the form. 



Teff{color)[WK] = J2 «i {colorf. (4.S) 

The fits were done for each set of models separately and they are shown in Figure 
14.81 with solid lines. The coefficients of the fits, ranges of validity for {B — V) and 
{V — I) colors, and dispersions are given in Table S]2j The fits to the E96 models are 
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Figure 4.8: Effective temperature versus (B — V) (left) and (V^— /) (right) from the 
models of E96 (filled circles) and DOS (open circles). The solid lines show polynomial 
fits for each set of models. 

characterized by a scatter of ~500 K in [B — V) and ~350 K in (\^ — /). For the DOS 
models the scatter is ~670 K in (i? — V) and ~800 K in {V — I). We do not have 
any strong argument to rule out either set of models. We therefore keep both results 
even if they show significant systematic differences. But we notice that the value 
of Te// during the recombination phase (where it is nearly constant) appears to be 
somewhat underestimated {Teff ~ 4600 K) by E96. We recall that these calibrations 
are only valid until the end of the plateau phase. 

In Figure 1^19) we show how these calibrations work for estimating Teff for SN 
1987A, SN 1999em, and SN 2003hn. As a comparison we included in these plots 
color temperatures obtained from the BB fits described in section 14.1.11 Note that 
for the three SNe the color temperatures are greater than T^ff, which is expected for 
"dilute" atmospheres whose continuum opacity is dominated by electron scattering. 
As said above, the usefulness of these fits is that they allow one to obtain Teff for 
any SN Il-P from their (B — V) or {V — I) colors, and thereby use this quantity to 
compare with hydrodynamical models. Note that we could equivalently have chosen 
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Table 4.2: Coefficients of the fits to Tg// (color) 







B-V 


B-V 


V-I 


V-I 


Hi 




E9Q 


D05 


E9Q 


D05 


ao 




0.790 


0.884 


0.957 


1.106 


ai 




-1.856 


-2.340 


-2.254 


-1.736 


a2 




4.055 


6.628 


5.922 


-6.403 


as 




-3.922 


-8.456 


-18.476 


33.762 


a4 




1.368 


4.619 


36.058 


-48.260 


as 






-0.849 


-25.291 


22.362 


rang 


;es 


[-0.2,1.15] 


[-0.2,1.15] 


[-0.1,0.65] 


[-0.07,0.83] 


rms 


[K] 


500 


670 


350 


800 



nff{color)[W^K] = ^:Lo«. (colorY 



to calibrate Rph vs color. 



4.2 Photospheric velocities 



Another important parameter to compare with observations is the photospheric ve- 
locity yielded by our models. Expansion velocity estimates from the minimum of 
several spec t ral li nes (Hq,, Fe ll A5169, and H^) of our data set were given by 
(120091). Since each line forms at a different shells, it is not straight for- 



Jones et al 



Jones et al 



ward to compare the observed velocities with photospheric velocities 
(120091 ) got around this problem and derived congruent calibrations between the veloc 
ity derived from the absorption minimum of such lines and the photosph eric velocity 



Jones et al 



using two independent atmosphere models (E96 and D05). As noted by 
(120091 ) provides a very good proxy to the photosphere velocity as it is not highly 
saturated as H^, and is present over most of the SN evolution. Ther efore, in the 
compa risons with our models, we chose to use the calibration given by iJones et al. 



(120091 ) for this particular line. However, in spite of the satisfactory behavior of this 
calibration in the high- velocity regime {v > 5000 km s~^), for lower velocities the 
atmosphere models fail to reproduce the behavior shown by the data (see Figure 9 of 



Jones et al. 



20091 ). Thus, we decided to complement our comparison with velocities 
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Figure 4.9: Effective temperatures for SN 1987A (upper panel), SN 1999em (mid- 
dle panel) and SN 2003hn (bottom panel ) as a function of time since explosion 
calculated using the polynomial fits given in table 14.21 for the models of E96 (open 
circles) and DOS (filled circles). As a comparison we include in these plots the color 
temperatures obtained from the BB fits described in section I4.1.1I 
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estimated from the Fe ll A5169 line for the case where it is present at later phases 
of the SN evolution and therefore the velocities are below the limit imposed by the 
calibration of hydrogen lines. 

For consistency with the atmosphere models, we define the photospheric posi- 
tion for our hydrodynamical models as the layer where the total continuum optical 
depth is r = 2/3. Note that to calculate r we do not include the opacity floor (see 
discussion in section I3.1.2p because such minimum is only included to take into ac- 
count line effects, and therefore is a bound-bound opacity that does not contribute 
to form the continuum. With this definition, the photosphere follows the recombi- 
nation wave, as we show in section 15.1.31 due to the fact that electron scattering 
is the dominant source of opacity (see Figure 13. ip . Had we included the opacity 
floor in the definition of r, the r = 2/3 surface would be pushed well above the 
recombination front. 



We remark that with this definition, the photosphere is essentially the surface 
of last scattering. However, the surface where the continuum is actually formed is 
located in a deeper layer called the "thermalization depth" and thi s is due to the 
dominance o: 



Hoflich 



1991 



the electron scattering oy er the absorption processes ( ISobolev 



1980; 



Montes &: Wagonerl . Il995[ ). With our simple prescription of radiative 



transfer we cannot accurately determine the location of the thermalization depth 
because there radiation decouples from matter. Note that the color temperature, 
determined by a black body fit to the broad-band photometry, is nearly coincident 
with the temperature of the thermalization depth but is greater than the effective 
temperature (as also shown in Figure 149]) . 
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Chapter 5 

Hydro-Model for SN 1999em 



In this Chapter, I present a detailed study of the prototype SN 1999em using the 
hydro dynamical code described in Chapter [3l The choice of this object was moti- 
vated by the fact that it is one of the best observed SNe of its type, both in terms of 
wavelength coverage and temporal sampling. In § 15.11 1 present the bolometric light 
curve and photospheric velocities for this object resulting from model calculations. 
A remarkably good agreement with observations was obtained using the following 
physical parameters: E = 1.25 foe, M = 19 M©, R = 800 Rq and = 0.056 Mq. 
In this analysis, 1 find that an extensive mixing of ^^Ni is required in order to re- 
produce a plateau as fiat as that shown by the observations. In § 15.21 1 discuss 
how our results compare with two pre vious hydrodynamical studies of SN 1999em 
( Baklanov et al. , 2005; Utrobinl . 2007 ). Finally, in § 15.31 1 study the possibility of 



fitting the observations with lower values of the initial mass consistently with the 
upper limits that have been inferred from pre-supernova imaging of SN 1999em in 
connection with stellar evolution models. 1 cannot find a set of physical parameters 
that reproduce well the observations for models with pre-supernova mass of < 12 
Mq, although models with 14 Mq cannot be fully discarded. 
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5.1 High- Mass Model of SN 1999em 



SN 1999em was dis covere d shortly after explosion o n 1999 October 2 9 UT by the 



Jones et al 



we assume 



LOSS program (JLj, Il999f ). Based on EPM study by 
SN 1999em exploded three days before discovery. This is consistent with the con- 
straint imposed by a negative detection (limiting magnitude 19) on an image obtained 
on 1999 October 20 UT. In order to compare with our model, we correct the times 
elapsed since explosion by time dilation based on the redshift of the host ga laxy. We 

(bood) to com- 



Leonard et al 



adopted the Cepheid distance of 11.7 Mpc given by 
pute bolometric luminosities as explained in § l4.1.2l using photometric data obtained 
at Cerro Tololo Inter- American Observatory (CTIO), Las Campa nas Obs e rvato ry 
(LCO), and the European Southern Observatory (ESO) at La Silla ( jHamuyl . 1200 ll ). 



We use a pre-supernova r nodel with init ial parameters consistent with the opti- 
mal hydrodynamical model of lUtrobinI ( 120071 ). Specifically, we adopt an initial mass 
of 19 Mq, radius of 800 Rq, and explosion energy of E = 1.25 foe (1 foe= 1 x 10^^ 
erg). This energy was released as thermal energy near the core of the object in a very 
short time scale as compared with the hydrodynamic time scale of our model. We 
also assume a nickel mass of 0.056 Mq which was determined by the luminosity of 
the radi oactive tail. Th is parameter is quite different to the nickel mass of 0.036 Mq 
used by lutrobinI f 20071) whi ch wa s determined using the quasi-bolometric luminosity 
given by Elmhamdi et al.l (2003). Their luminosity was based on the integration 
of only UBVRI photometry with a constant value of 0.19 dex added to take into 
account the infrared luminosity, while our bolometric luminosity is based on a quan- 
titative study of the bolometric correction for SNe II-P (see §|4]). No te also that our 



Baklanov et al. 



value for the nickel mass is closer to the estimate of 0.06 Mq given by 
(120051 ). In the following analysis, we denote this model as m19r8e1.25ni56. 



In our calculations we remove the central 1.4 Mq which is assumed to form a 
neutron star. The initial density profile as function of mass and radius for model 
m19r8e1.25ni56 is shown in Figure ISTTl Note that the initial structure is composed 
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of a dense core and an extended envelope characteristic of a red supergiant. The 
chemical composition profiles are shown in Figure 15.21 The outer parts of the enve- 
lope, M > 9 Mq, have a homogeneous composition with mass fractions for hydrogen 
of X = 0.735, for helium of F = 0.251, and a metallicity of Z = 0.02. From there in- 
ward, hydrogen and helium are mixed in order to prevent a sharp boundary between 
the H-rich and the He-rich layers. Such sharp boundaries are characteristic of stellar 
evolution models but fail to reproduce the observations. Note that we allow H to 
mix very deep inside the core and ^^Ni is mixed out in the envelope up to ~15 Mq. 
This type of mixing is presumably due to Rayleigh- Taylor instabilities that occur 
behin d the shock front, as obtained in mu lti-dimensional hydrodynamic calcula- 



tions (jMueller et 



(IShigeyama et al 



199 



1988 



Kane et al 



Woosley et al. 



2000h and supported bv studies of SN 1987A 



1988 



Arnett 



1988 



Blinnikov et al. 



12000 



among others). The presence of H in the core leads to a smooth transition between 
the plateau and the radioactive tail. The distribution of ^^Ni to external layers helps 
to reproduce a plateau as flat as that observed in SN 1999em (see section 15.1. 3p . 




Figure 5.1: Initial density distribution with respect to interior mass (left) and radius 
(right) for the pre-supernova model m19r8e1.25ni56. 



Figure 15.31 shows a comparison between the bolometric light curve (solid line) 
obtained with our code for model m19r8e1.25ni56, and the observations of SN 1999em 
(dots). The luminosity due to ^^Ni — ^^Co — ^^Fe (dashed line) is also shown. Note 
the very good agreement between model and observations. The largest difference ap- 
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Figure 5.2: Abundance distribution in the pre-supernova model m19r8e1.25ni56 
with respect to Lagrangian mass. For clarity, the abundance of carbon (C) and 
oxygen (O) were multiplied by 3 and the abundance of ^^Ni, by 20. Note that the 
^^Ni is uniformly mixed in the outer envelope until 15 Mq. 



pears during the earliest phase and the transition to the radioactive tail. At earliest 
epochs our bolometric corrections have the highest uncertainties, as noted in § I4.1.2[ 
and during the transition between the plateau and the radioactive tail, the diffusion 
approximation breaks down because the object is almost completely recombined and 
the photosphere is not well defined. During the tail, the bolometric luminosity is 
completely determined by the luminosity of radioactive decay, which makes our cal- 
culations more reliable. Although not shown here, our code shows that the shape of 
the light curve at the end of the plateau is very sensitive to the properties of the core, 
such as mixing and the form of the density transition between the helium-rich core 
and the hydrogen-rich envelope. In principle, it is possible to find an even better fit 
to the observations but, given the uncertainties introduced by our approximations, 
such a detailed study should be considered as meaningless. 



Figure 15.41 shows the photospheric velocity evolution of our model compared 
with observed photospheric velocities as explained in § 14. 2[ We have also included 
in the plot the spect roscopic velocities measured from the absorption minimum of 
the Fe ll A5169 line (jJones et al.l . |2009| ). reveling that the velocities estimated from 
this line samples quite well the photospheric velocity at all epochs when it is present. 
However, it is not present at the very earliest times {t < 10 days). Note the very good 
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Figure 5.3: The bolometric light curve for model m1 9r8e1.25ni56 (sol i d line ) com- 
pared with the data of SN 1999em as calculated by iBersten fc Hamuyl (120091 ) (blue 
dots). The luminosity due to the ^^Ni — )■ ^^Co — )■ ^^Fe decay is also shown (dashed 
line). The physical parameters used in the model are indicated. 



overall agreement between model and observations. Fe ll velocities match quite well 
the model photospheric velocities except at the latest times. However, it is important 
to remark that the photospheric velocity is to be considered a good discriminator 
between models only at the early phases of the evolution while the object is not 
completely recombined. This is because, line velocities become poor photospheric 
velocity indicators with time, and the photosphere begins to lose its meaning in 
our models as the ejecta becomes nearly completely recombined at the end of the 
plateau. 

As shown in Figures 15.31 and 15.41 we obtain a remarkably good agreement 
with observations (bolometric light curve and photospheric velocity evolution) for 
SN 1999em despite the simplifications used in our code. These results are very 
encouraging and give us confidence in the capability of our code to infer physical 
parameters and to study their effect on the observed quantities with the ultimate 
aim of understanding the physics of SNe II-P. 

In the following subsections, we describe in some detail the evolution of the 
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Figure 5.4: Evolution of the expansion velocity of the photosphere for model 
m19r8e1. 25ni56 (soli d line) compared with observed photospheric velocities cal- 
culated by I Jones et al.l (120091 ) using a calibration between the velocity derived from 
the absorption minimum of H/j and two atmosphere models: E96 (filled squares) and 
DOS (open circles). We also include the velocities calculated from Fe ll A5169 (filled 
circles). 



SN for this particular model. At a glance, the LC is distinguished by three phases: 
a) an outburst followed by strong cooling, b) a plateau, and c) a radioactive tail. 
Each phase is essentially determined by the interplay between the main heating and 
cooling mechanisms. We thus focus our discussion on these processes along the SN 
evolution. As a reference. Table ETT] gives a summary of some properties of the model 
at characteristic times during the evolution. The initial phase due to shock breakout 
is discussed in § I5.1.H the adiabatic-cooling phase where the homologous expansion 
is reached is addressed in § 15.1.21 the cooling and recombination phase is described in 
§ I5.1.3[ and finally the radioactive-decay processes that power the tail are explained 
in §[5X1 



5.1.1 Shock wave propagation and the early evolution 

A powerful shock wave (SW) begins to propagate outward through the envelope 
when we artificially inject energy near the center of the star (assumed to occur 
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Table 5.1: Properties of model m19r8e1.25ni56 at selected time of their evolution. 



Phase 


t 

[days] 


log{Lboi) 

[erg s^i] 


log{T,ff) 
[K] 


log(Tph) 
[K] 


[lO^^cm] 


[lO^cm s^l ] 


(%) 


(%) 


Peak 


1.5 


44.5 


5.05 


5.03 


0.54 


3.36 


0.04 


72 


Adiabatic Cooling 


5 


42.4 


4.10 


4.05 


3.8 


11.5 


0.2 


93 


Plateau 


50 


42.1 


3.73 


3.77 


14.7 


3.47 


0.65 


99 


Transition 


114 


41.7 


3.74 


3.74 


10 


1 


1.2 


99.7 



^ Erad is defined by Jq* Ljoj dt and En is the total kinetic energy in the mass motions. In the table we show 
these quantities normalized to the initial injected energy of model m19h8e1.25ni56, Eq = 1.25 foe. 



at t = 0). This energy is initially released as internal energy, and part of it is 
rapidly transformed into kinetic energy (e.g., by t = 0.5 days, the total energy is 
approximately equally divided between kinetic and internal energy). The velocities 
acquired by matter are so high that they exceed the local speed of sound, leading 
to the formation of a SW. The SW heats and accelerates the matter depositing 
mechanical and thermal energy into successive layers of the envelope until it reaching 
the surface, where photon diffusion dominates the energy transfer, and energy begins 
to be radiated away. 

Figure 15.51 shows the effect of the SW propagation on different physical quanti- 
ties (velocity, density and temperature) inside the star. The shock front is evidenced 
as a sudden change in these quantities. As the shock moves outward, the material 
behind is accelerated and heated. Note that the outermost layers, with a sharp de- 
cline in density, acquires very high velocities. It represents a small fraction of the 
star. It is also clear from the velocity profiles that the inner parts of the object are 
decelerated. This deceleration is due to the interaction of the dense core with the 
extended hydrogen-rich envelope. The bottom right figure shows changes in radius 
for different shells. From this we can deduce the location of the shock front at any 
time as that of the innermost shell which has constant radius. 

At t = 1.36 days, the SW reaches the surface of the object, which produces the 
first electromagnetic manifestation of the explosion (although neutrinos and gravita- 
tional waves escape well before). The effective temperature and bolometric luminos- 
ity suddenly rise and reach their maximum values a few hours after breakout, specif- 
ically, at t = 1.47 days with values of Lpeak = 3 x 10"^^ erg s~^ and Tpeak = 1.1 x 10^ K 
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Figure 5.5: Changes in density (top left), temperature (top right), and velocity 
(bottom left) profiles as a function of interior mass during shock propagation for 
model m19r8e1.25ni56. Some of the curves are labeled with the time elapsed 
since the energy is injected. The initial density profile is also shown [t = 0; blue 
dashed line). Note that a very small amount of material near the surface is strongly 
accelerated as the SW passes through the steep density gradients present in the 
outermost layers. Bottom right: temporal evolution of the radial coordinates 
corresponding to different interior layers. 
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as shown in Figure 15.61 At these temperatures, the peak of the emitted spectrum is 
in the UV range. 



Hereafter, the star begins to expand and cool very quickly, leading to an increase 
in photospheric radius and a decrease in temperature in the external layers. The 
bolo metric luminosity abruptly decreases but, according to the decrease in effective 
temperature and consequent shift of the emission peak to longer wavelengths, the 
luminosity in the optical range increases. As a result, a sharp peak in bolo metric 
luminosity and temperature is produced, as shown in Figure 15.61 In our model 
we obtain a decrease of 1.5 dex in luminosity only 6.8 hours after peak brightness. 
During this time the total energy radiated is 2.2 x 10^^ erg, emitted essentially as 
a UV flash. The short duration of the breakout explains why so few SNe II-P have 
been observed during this phase. 
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Figure 5.6: Bolometric luminosity (left) and effective temperature (right) during 
shock breakout 



At temperatures as high as those left by the passage of the shock wave, the 
stellar matter is completely ionized which implies that the breakout is accompanied 
by a strong increase in opacity. The position of the photosphere during the outburst 
nearly coincides with the outermost shell. This behavior continues until the onset of 
recombination. Therefore, at early stages previous to recombination, the velocity of 
matter in the photospheric position samples the very high velocities of the outermost 
layers and reaches values close to 1.2 x 10^ km s~^. 
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From the energetics point of view, by t = 0.5 day the total energy is ap- 
proximately divided in equal proportions between internal, dominated by radiative 
contribution, and kinetic components. Short after breakout the kinetic energy com- 
pletely dominates the energetics (the energy radiated away at any given time is less 
than 1.5% of the total energy). 

It is possible to estimate the average velocity during the SW propagation: given 
that the SW takes 1.36 days to emerge and considering a radius of i? = 5.5 x 10^^ cm 
(800 Rq), we obtain an average speed of v = 4680 km s~^. It is also interesting 
to calculate the e xpect ed time for breakout using the analytic expression given by 
Shigeyama et al. 



thk ^ 1-6 



Rn 



5Oi?0 



IOMq 



E 



1 X 10^^ erg 



1/2 



hr. 



(5.1) 



where Rq is the initial radius, Mej is the ejected mass, and E is the explosion energy. 
Using the values for our initial model we obtain tbk = 1-26 days, in very good 
agreement with our numerical calculation. 



5.1.2 Adiabatic cooling and homologous expansion 

The breakout is followed by a violent expansion, resulting in the cooling of the 
outermost layers. During expansion, only a small fraction of the photon energy can 
diffuse into the surroundings. Therefore, it is possible to consider the cooling process 
to be approximately adiabatic and this approximation remains valid while the time- 
scale for radiation diffusion is much longer than the expansion time-scalqj. Note that 
there are two mechanisms to cool a SN: (a) loss of photons or diffusion cooling, and 
(b) its own expansion or "adiabatic cooling". If one of these processes dominates 
then we say that the cooling is carried out by such process. In our model, more than 



^The expansion time-scale, Th = R/v, increases with time while the diffusion time-scale, tj — 
KpR^/c, decreases because p oc R"^. Thus, there is a time after which the condition for this 
approximation breaks down. 
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90% of the decrease in internal energy happens as adiabatic coohng during the first 
18 days after explosion, when the first layers with neutral hydrogen appear. After 
that, diffusion cooling begins to be significant and the adiabatic cooling phase comes 
to an end. 

The internal temperature (and internal energy) decreases almost adiabatically, 
i.e. proportional to due to the dominance of the radiative term, and quickly 
reaches a value near the recombination temperature of hydrogen. Quantitatively, the 
effective temperature goes from values close to 10^ K at the time of the burst when 
the matter is totally ionized, to T^s = 10'* K five days later. At these temperatures 
hydrogen begins to recombine. At the same time, the luminosity reaches L = 1.9 x 
]^q42 gj,g g-i photospheric radius rapidly increases to i?ph = 5.2 x 10^^ cm. 

After that, the luminosity decreases slightly as a result of the slower decrease in 
temperature and continuous increase in radius. By day 18, when the first layers 
of neutral hydrogen appear, we have T^s = 6960 K, L = 1.2 x 10"^^ erg s~\ and 
Rph = 9.1 X 10^^ cm. Note the slight change in luminosity, of only 0.14 dex, between 
day 7 and 18 (Figure 15731) . 

A few days after breakout, the acceleration of the material comes to an end and 
the expansion becomes homologous. The matter reaches velocities of the order of 
1.2 X 10^ km s~^ in the outermost layers and the object enters a state of free expansion 
where forces of pressure and gravitation do not have any dynamical effect on the 
system. The homologous regime is characterized by a constant velocity in each layer, 
a linear growth of the radial coordinate with time (r oc t) and a density distribution 
decreasing with time as p oc t~^. This behavior is clearly shown in Figure [STSl for t > 3 
days and in Figure 15.71 The nearly constant shape of the density and temperature 
profiles at late times are a result of the expansion being approximately homologous. 
This is also evident in the linear behavior of the radial coordinates for different 
mass shells. The condition of constant velocity is very nearly satisfied for each shell, 
although the transition between acceleration and homologous expansion happens at 



slightly different times for different mass shells (see Figure I5.7ji n. 

^After day 8, the photosphere remains located in shells where homology has been reached. 
Therefore, our models show that the expanding photosphere method (EPM) which is used to 
estimate distances assuming homologous expansion, can be safely applied at epochs later than this. 
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Figure 5.7: Evolution of the velocity for different layers (solid lines) and the velocity 
of material at the photospheric position (dashed line). Note that the asymptotic 
constant velocity is reached at latter times for the more external shells. 



As mentioned in § I5.1.H at the time of breakout the internal energy of the 
envelope is small compared to the kinetic energy. Moreover, the object expands by a 
factor of ~ 17 during this phase to a radius of i?ph ~ 9 x 10^^ cm, before it becomes 
transparent. Thus, most of the internal energy left by the passage of the shock has 
been already degraded by the adiabatic expansion, and only a small residual will be 
released as observable radiation in the following phase. 



5.1.3 Cooling and recombination wave 



The appearance of regions with neutral hydrogen determines the onset of a recombi- 
nation wave (RW). The duration of the RW is related to the total hydrogen mass and 
how deep hydrogen has been mixed into to the progenitor object. As time goes on, 
hydrogen recombination occurs at different layers of the object as a wave propagates 
inwards (in Lagrangian coordinate; see Figure [5l8| top). Because the opacity is dom- 
inated by electron scattering, it decreases abruptly outwards at the recombination 
front (Figure ESI top right), increasing the transparency of these layers and allowing 
the radiation to easily go away. Consequently, the internal energy (mostly of the 
radiation field) is efficiently radiated away, and the temperature drops sharply from 
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~10000 K to ~5500 K at the recombination front (Figure [5.8^ bottom). In other 
words, a coohng wave associated with the transparency and induced by a recombi- 
nation wave propagates inward through the envelope. This is usually called "cooling 
and recombination wave" (CRW). 

In our model the recombination begins approximately at day 18. From them 
on, a CRW develops which moves inward in mass until all the matter is completely 
recombined by day 130 when the recombination front reaches the innermost H-rich 
layers of the envelope (at m =2 Mq in our model), and the luminosity suddenly 
drops, defining the end of this phase. The propagation of the CRW is clear in 
Figure 15.81 where the evolution of the fraction of ionized hydrogen and temperature 
profiles as a function of mass for selected times are shown. Note also the behavior of 
the opacity (see Figure [578| top right) which depends strongly on the ionization state 
of the matter. In our model, a strong drop in the opacity of the outer layers is seen at 
day ~18, which leads to a considerable decrease in the optical depth of these layers 
and causes an inward drift in mass of the photosphere — defined at a fixed optical 
depth of T = 2/3. The photosphere begins to follow the CRW, as shown in Figure 
15.81 It is important to note that we define the position of the photosphere using only 
continuum opacity sources, i.e. excluding the opacity floor (see section . 

The photosphere, as defined here, is nearly coincident with the outer edge of 
the CRW. Therefore it is just the location of the photosphere with respect to the 
CRW what sets the value of the photospheric temperature close to the temperature 
of hydrogen recombination (T ~ 5500 K). However, note that the effective temper- 
ature does not necessarily have a constant value. Since the effective temperature 
is defined as Tjj = L/ttctR'^^ and the luminosity remains nearly constant outside 
the recombination front (see Figure IS.lOp . any changes in effective temperature are 
related to changes in photospheric radius. 

The CRW divides the object in two distinct regions: (a) an inner zone which is 
hot, optically thick and ionized, and (b) an outer zone which is relatively cold, opti- 
cally thin and completely recombined. Matter in the inner zone is opaque: radiative 
transfer is too inefficient to produce any appreciable flow of radiative energy (this 
would be strictly fulfilled if ^^Ni were confined to the innermost layers; see the dis- 
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Figure 5.8: Evolution of the fraction of ionized hydrogen (top left), opacity (top 
right), and temperature (bottom left) as a function of mass. The time elapsed 
since the energy is injected and the photospheric position (blue dot) are indicated for 
each curve. The photosphere is nearly coincident with the outer edge of the CRW 
which moves inward in mass coordinate. Note also that the opacity in the outer 
shells is nearly constant with a value of 0.4 cm^ g~^, close to the electron-scattering 
opacity for matter composed of pure hydrogen, above which it suddenly drops. 
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cussion below) and the matter cools down almost adiabatically. On the other hand, 
the external layers are transparent and practically do not radiate. Therefore, it is 
within the CRW interphase where almost the entire radiant flux is released. When 
matter passes through the CRW, particles begin to be cooled by radiation, emitting 
more light than they absorb, and the radiant flux increases. This way, the radiative 
flux emerging from the CRW front carries away internal energy of the matter that 
is cooled by the wave . More details on the properties of the CRW are given by 



Grassberg et al.l ( 1l97ll ). 



Note that the bulk of the radiation does not diffuse to the photosphere; the 
photosphere instead moves inward, allowing the radiation to escape sooner than it 
would for a photosphere fixed at the outer boundary of the ejected mass. That is, 
the recombination process is responsible for the energy release during this phase. 
It should be noted that most of this energy comes from the energy de- 
posited by the SW and not from the recombination itself. In order to test 
the previous statement, we ran our code without including the energy released by 
recombination of ions with electrons. The result was that no appreciable change ap- 
peared; quantitatively, the differences between both calculations are less than 0.04 
dex during all the evolution (see Figure [5^ . 

Note, however, that in model ]Vi19r8e1.25ni56 where we assumed an extended 
^^Ni mixing, the flux of energy inside the CRW is not negligible (see left panel 
of Figure IS.lOp . This is due to the fact that the photosphere meets regions with 
radioactive material earlier than in the case where ^^Ni is confined to the innermost 
layers. Therefore, the energy deposited by radioactive decay provides additional 
power to the LC. A different behavior is seen when ^^Ni is confined to the innermost 
layers (inside 2.5 M©; see right panel of Figure IS.lOp . In this case, the statement 
that the energy flux inside the CRW is very small is fulfilled at least until day ~120 
when the radiative diffusion of the radioactive decay from the central layers begins to 
dominate. At earlier times, it is possible to see the outward diffusion of radioactive 
energy. Note also that for both models the luminosity outside the front is nearly 
constant. 

The assumption of extended mixing was necessary in order to obtain a plateau 
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Figure 5.9: Comparison between bolometric light curves with (sohd hne) and without 
(dashed hne) energy released by recombination. It is clear that the recombination of 
ions with electrons does not contribute significantly to the light curve. The differ- 
ences between both calculations are less than 0.04 dex during whole the evolution. 



as flat as the one observed for SN 1999em. Figure 15.111 shows a comparison of 
the bolometric LC for three cases: 1) with extended ^^Ni mixing, 2) no mixing, 
and 3) without ^^Ni. For case 1), Ni begins to affect the LC by day ~35 while 
for case 2) the effect of nickel heating is delayed until day ~75. On the other 
hand, in the former case there is a less direct energy deposition at the center of 
the object, and the plateau declines earlier and steeper to the tail than in case 2). 
Our assumption of mixing of ^^ Ni into the hydrogen envelope is not unreasonable as 



show n in studies of SN 1987A (IShigeyama et al 



1988; 



1988; 



Blinnikov et al 



Woosley et al. 



1988 



Arnett 



2000l among others). This type of mixing is presumably 



due to Rayleigh- Taylor instabilities that occur be hind the shock front , as obtained 



1991 



Kane et al.l . 



in ni ulti-dimensional hydrodynamic calculation s (|Mueller et al.l . 
2000). However, it is important to note that UtrobirJ ~ 2007 ) found an excellent 



agreement with observations of SN 1999em by conflning ^^Ni to the innermost layers 
and good agree r nent with observations of SN 1987A assuming moderate ^^Ni mixing 
Jutrobinl . ll993l . l2004l ). 
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Figure 5.10: Evolution of the interior luminosity as a function of mass. The labels 
indicate the time since explosion. Left: ^^Ni is mixed into the hydrogen-rich envelope 
up to 15 Mq. Right: ^^Ni is confined to the layers inside 3.5 Mq. The inward 
propagation of the recombination front is clear in both cases. In the mixed case, 
the recombination front propagates slower than in the unmixed one because the 
temperature near the front is higher due to radioactive heating. Note that for the 
unmixed case the outward diffusion of the radioactive energy is clear while for the 
mixed case this effect occurs too early to be noticeable. 
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Figure 5.11: Comparison between bolometric light curves for extended ^^Ni mixing 
(solid line) and a model with mixing of ^^Ni up to 3.5 Mq (short-dashed line). We 
also include the case of a model without ^^Ni (long-dashed hne) in which case the LC 
falls abruptly after the plateau phase. The presence of ^^Ni extends the plateau and 
increases the luminosity. This is essentially produced when the CRW reaches layers 
with ^^Ni which can thereby power the LC directly. The extended ^^Ni mixing reveals 
this effect earlier in the evolution, which produces a fiat plateau in concordance with 
the observations of SN 1999em. 
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The evolution of the velocity of matter at the photospheric position fph, photo- 
spheric radius (-Rph) and mass above the photosphere (Mph) are shown in Figure [5rT2l 
Initially, the photospheric velocity evolves rapidly, the photospheric position follows 
a linear behavior, and there is very little mass above the photosphere. Later on, 
when the recombination sets in at t ~ 18 day, Rp^ begins to differ noticeably from 
the linear behavior, the mass above the photosphere increases, and fph decreases 
because it samples increasingly inner, slower material. Finally, when all the matter 
is recombined at t ~ 110 days, -Rph and fph sharply turn down and the luminosity 
undergoes a rapid decrease. Note, however, that the photospheric radius does not 
drop immediately to zero. This is due to the heating caused by radioactive decay 
which produces some ionization of the gas. On the other hand, the luminosity un- 
dergoes a rapid decrease to values close to the luminosity of radioactive decays. If 
there was no ^^Ni (case 3) the SN luminosity would abruptly vanish at this point, as 
shown with a long-dashed line in Figure 15.111 




Figure 5.12: Evolution of photospheric radius (Rph) in units of 10^^ cm, photospheric 
velocity (V^h) in units of 5 x 10^ cm s~^, and mass above the photosphere (Mph) in 
units of Mo (Mo = 19 Mq) for model m19r8e1.25ni56. 



In conclusion, the CRW has a duration of ~100 days. The luminosity experi- 
ences small changes during the CRW propagation. In order to produce a plateau as 
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flat as that observed in SN 1999em, we need to invoke significant Ni mixing in the 
H-rich envelope. Thus, the plateau can be seen as a combination of CRW 
properties plus some additional energy provided by radioactivity. Defin- 
ing, from an empirical point of view, the plateau phase as the period of time when 
the luminosity remains constant within 0.5 mag of the value at day 50, its duration 
is 103.5 days, and the total energy emitted during this phase is ~1.12 x 10^^ erg. 
Note, that with such definition the plateau phase includes the final stages of the 
early adiabatic cooling phase and not only the recombination phase. 

Radioactive tail 

The late behavior of the LC (at t > 130 days) is dominated by the energy released 
from radioactive decay. Without radioactive material, the luminosity would abruptly 
vanish when hydrogen gets completely recombined as shown in Figure 15.111 Instead, 
the observed LC decreases to values close to the instantaneous rate of energy depo- 
sition by the radioactive decay, 

L ~ 1.43 X IO^^Mni/Mq exp(-t/111.3). (5.2) 

The decline of the LC at t > 130 days is nearly linear in concordance to the expo- 
nential decline of the radioactive decay law. That is, the diffusion time for optical 
photons becomes so small that the decay energy is radiated away nearly instanta- 
neously, while the gamma-ray optical depth is still sufficiently large in order to allow 
a nearly complete local deposition of the decay energy. The bolometric luminosity in 
this part of the LC is a direct measure of the Ni mass synthesized in the explosion. 
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5.2 Comparison with other hydro- models of SN 
1999em 



In this Section, we compare our results with those of previous hydrodynamical stud- 
ies of SN 1999 e m. T here are two previous such studi es of SN 1999e m: one by 



Baklanov et aL 



fl2005r ) (BBP05 hereafter) and another bv lUtrobin] fl2007f ) (U07 here- 
after). In Table 1521 we summarize the physical parameters obtained for SN 1999em 
in the three works along with the distance, explosion time, and the chemical com- 
position assumed in each model. It is interesting to note that the three studies 
employ initial density distributions and chemical composition distributions of non- 
evolutionary models. In our work we assume H and He profiles very close to those 
adopted by U07, yet a very different ^^Ni distribution. While we assume an extended, 
uniform mixing of Ni (see Figure |5T2|) . U07 confined ^^Ni to the innermost layers (see 
their Figure 2). In the case of BBP05, H and He were assumed to be uniformly 
mixed throughout the envelope and a radial distribution of ^®Ni up to ~15 Mq was 
adopted (see their Figure 9). Therefore, our ^^Ni distribution and the resulting ^^Ni 
mass are closer with those of BBP05 than those of U07. We also note that the mass 
of the compact remnant — which is left aside from the calculations — is different in 
all three works. We assume a compact remnant of 1.4 Mq while U07 assumed 1.58 
Mq and BBP05 removed everything within a radius of Rc = 0.1 i?©. 



Table 5.2: Comparison between physical parameters for SN 1999em from three dif- 
ferent hydrodynamical codes 





D 


to 






E 


M 


R 


Mm 


Ni mixing 


Code 


[Mpc] 


[JD-2451000] 


^sup 


Z 


[foe] 


[Mq] 


[Rq] 


[Mq\ 


[Mq] 


This work 


11.7 


477.90 


0.735 


0.02 


1.25 


19 


800 


0.056 


~15 


BBP05'' 


12 


468.90 


0.7 


0.004 


1 


18 


1000 


0.06 


-15 


U07'' 


11.7 


476.90 


0.735 


0.017 


1.3± 0.1 


20.58±1.2 


500± 200 


0.036±0.009 


~2.5 



^ iBaklanov et al.l | |2005|) 



^ lUtrobid ||2007|1 



The approaches used in each work are very different. U07 used a hydrody- 
namical code with a one-group approximation for the radiative transport, including 
non-LTE treatment on opacities and thermal emissivity, non thermal ionization. 
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and expansion opacity. Their calculations, although involving a more sophisticated 
radiative transfer treatment, yielded bolometric light curves comparable to those 
given by our code. However, the procedure to produce the observed bolometric light 
curve for SN 1999em was quite different in both cases. While U07 based their cal- 
culations on the integration of UBVRI photometry with a constant value of 0.19 
dex added to take into account the infrared luminosity, we derived UBVRIJKLM 
bolometric lumino s ities f rom color- dependent bolometric corrections. As noted by 
Bersten fc Hamuvl koO^ . the inclusion of the L and M bands in the calculation 



of the bolometric correction during the tail phase is very important and it most 
likely leads to the difference in the nickel mass estimated in this work and in U07. 
BBP05, in turn, used a multi-group hydrodynamical code which allows to calculate 
LC in different photometric bands. This code also included the expansion opacity 
effect. Therefore, BBP05 were able to compare their model with UBVRI LC sep- 
arately without being affected by the uncertainties in the calculation of bolometric 
luminosities. 



Table As shown Table l5^ in the parameters yielded by our calculations are 
intermediate between those estimated by BBP05 and U07. The largest differences 
are 2.6 Mq in mass, 500 Rq in radius, 0.3 foes in energy, and 0.045Mq in ^^Ni 
mass. It is quite satisfactory that our simple prescription for the radiative transfer 
yields physical parameters similar to those obtained from more sophisticated codes, 
considering that (a) the models and methods used in each study are quite different, 
(b) there are differences in assumed quantities, such as distance, explosion time, 
mass cut, photometry, and (c) there may be a degree of degeneracy among explosion 
energy, initial radius and initial mass within each model. 



5.3 Low Mass Model of SN 1999em 



The mass of the progenitors of SNe II-P can be derived from hydrodynamical model- 
ing of light curves and expansion velocities, or from the detection of the pre-supernova 
object in archival images of the host galaxy in connection with stellar evolution mod- 
els. At the moment three progenitor stars of SNe II-P have been firmly detected and 
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there have been negative detections for 17 other object s which have led to upper 
hmits of the progenitor star masses (jSmartt et al.l . 120091 ) . Among these, only three 
SNe II-P have masses derived using hydrodynamical models: SN 1999em, SN 2004et, 
and SN 2005cs. In all of these cases the masses estimated from the hydrodynamical 
modeling are systematically higher than the values derived fro m the other method . 
These discrepancies have been noted previously in the literature (lUtrobin fc Chugail . 



2008; 



Smartt 



20091 ). For the particular case of SN 1999em, the pre-supernova images 



give an upper limit of 15 Mq for the progenitor star in the ZAMS. It is expected that 
the mass of the pre-supernova object should be lower that this due to possible mass 
loss episodes during the evolution of the star. In this Section we use our hydrody- 
namical model to explore a low mass range, consistent with pre-supernova imaging 
of SN 1999em, in order to test how well we can reproduce the observed properties of 
this object. 



We calculate several models using two different values of the pre-supernova 
mass: M = 12M0 (see Figures EIIS and EUl]) , and M = 14Mq (see Figures EUS] 
and I5.16p . In all Figures, we show the bolometric light curves and photospheric 
velocities (fph) for different values of the injected energy (E), initial radius (R), and 
degree of ^^Ni mixing. Table |5T3] gives a summary of the parameters used. We remark 
that in all of these models, we have assumed the same value for the ^^Ni mass as 
that of model m19r8e1.25ni56, i.e. O.OSGM©, because this value is required to fit 
the tail of the LC. Also note that the degree of ^^Ni mixing is taken as a fraction of 
the initial mass of the object in order to compare the different degrees of mixing in 
a consistent way when models with different initial mass are used. For example, a 
model with a mixing of ^^Ni up to O.8M0 has an equivalent degree of ^^Ni mixing as 
that of model m19r8e1.25ni56. 



For the case of 12 Mq, we have calculated nine models. We first analyze 
the effect on the results of the variation of one parameter while keeping the other 
parameters fixed. In the upper panel of Figure 15.131 we show the effect on the LC 
and Vph of the variation of E for a model with R = 800 Rq and ^^Ni mixing up to 
0.8 Mq. In the bottom panels of Figure [5.13[ we show the effect of varing R, for a 
model with E = 1 foe and for the same ^^Ni mixing. The sensitivity of the LC and 
Vph on the extent of ^^Ni mixing is shown in the upper panel of Figure 15.141 for a 
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model with E = 1 foe and R = 800 R 



An examination of these Figures yields the following conclusions: 1) higher 
injected energy produces higher luminosity and a shorter plateau, 2) larger initial 
radius produces higher luminosity and a longer plateau, 3) more extended ^^Ni mixing 
produces higher luminosity and a shorter plateau, 4) ^^Ni mixing and initial radius 
have a small effect on the Vph evolution as compared with that of the explosion 
energy, and 5) as expected, there is no effect of any of these parameters on the 
tail luminosity. Note that the sensitivity of the LC on -E, R, a n d ^^Ni niixing is in 



qualitative concordance with analytic studies by 
pr evious numeri cal st udies by 



by lYouneJ ( 120041 ) and 



Utrobin 



Arnett 



(119801) IPopovI (119931 ). and 



Litviiiova fc Nadezhin] ( 119831 . Il985[ l and more recently 
J2OO7I ). 



Although we have considered several values of E, R, and ^^Ni mixing, neither 
of the models with M = 12 Mq gives a good representation of the observations of 
SN 1999em. Note that we have calculated also two other models, m12r15e05ni56 
and m12r15e08ni56, adopting even higher values of the initial radius. These mod- 
els are shown in the bottom panel of Figure 15.141 Despite the good representation 
of the LC provided by model m12r15e05ni56, it fails to reproduce the observed 
photospheric velocities. We conclude that it is not possible to reproduce the ob- 
servations of SN 1999em using models with a pre-supernova mass of M = 12 M©. 
Clearly, masses lower than this value would also fail to match the observations. 



For the case of 14 Mq, we have calculated ten models for different combinations 
of R, E and ^^Ni mixing. Specifically, we used three different values of the initial 
radius: R = 800, 1000, and 1200 i?0, four values of explosion energy: E = 1.1, 1.0, 
0.9, and 0.8 foe, and three different degrees of ^^Ni mixing: 0.2, 0.5, and 0.8 Mq 
(see Table 15. 3p . Figures 15.151 and 15.161 show the results of these models. At first 
sight we find that these models agree better with the observations than the 12 Mq 
ones. Therefore the case of 14 Mq deserves a more detailed analysis. As described 
below, we chose the grid of parameters based on the known effects of each physical 
parameter on the light curves and velocities, and trying to reach the best possible 
agreement with the observations. As noted in Section 15.11 in order to assess the 
validity of models we will focus on how well they reproduce the observed plateau 
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Figure 5.13: Comparison between models and observations of SN 1999em for a low 
value of the pre-supernovae mass of M = 12 Mq and the same value of ^^Ni mass and 
mixing than for the reference model m19r8e1.25ni56. (Left panels): bolometric 
LCs. (Right panels): photospheric velocity evolution. (Upper panels): models 
with different explosion energies as indicated in the labels, and a fixed value of 
R = 800 Rq. (Lower panels): models with different initial radii as indicated in the 
labels, and a fixed value of E = 1 foe. 
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Figure 5.14: Comparison between models and observations of SN 1999em for a low 
value of the pre-supernovae mass of M = 12 Mq and the same value of ^^Ni mass than 
for the reference model m19r8e1.25ni56. (Left panels): bolometric LCs. (Right 
panels): photospheric velocity evolution. (Upper panels): models with different 
mixme; of ^^Ni, E = 1 foe and R = 800 Rq. Note that the de Efree of ^^Ni is indicated 
in each figure as a fraction of the initial mass of the model. (Lower panels) : models 
with an initial radius of R = 1500 Rq and two different values for the explosion 
energy as indicated (corresponding to models ml2rl5e08ni56 and ml2rl5e05ni56; 
see Table 1531) . Model ml2rl5e05ni56 reproduces very well the bolometric light curve 
of SN 1999em, while it underestimates the photospheric velocity evolution. 
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luminosity and length, and the early photospheric velocities, but we will disregard 
differences in the early LC before the plateau phase where the uncertainties in both 
models and data are the largest. 

Our reference model is shown with a solid line in the upper panels of Figure [5. 151 
and has the following parameters: E = 1 foe, R = 800 Rq, and 0.8 Mq of ^^Ni mixing. 
We note that this reference model has the same initial radius and ^^Ni mixing as 
model m19r8e1.25ni56 for 19 Mq. The energy has been reduced in order to com- 
pensate the effect on the luminosity of using a lower mass. As shown in Figure 15.151 
the reference model produces the correct plateau luminosity and photospheric veloc- 
ities although the plateau duration is too short as compared with the data. With 
the aim of remedying this situation while keeping the mass fixed we invoked lower 
energies and larger radii (dashed and dotted lines). While this served to improve 
the issue of the plateau length, the comparison with photospheric velocities at early 
times became poorer as expected due to the lower energies. We therefore consider 
these models to be unlikely. In a further attempt to find a good 14 Mq model, we 
decided to vary the mixing of ^^Ni while keeping the other parameters as in the ref- 
erence model. The results are shown in the lower panels of Figure 15.151 We now note 
that while a reduction of the ^^Ni mixing serves to increase the length of the plateau, 
the shape of the LC and the plateau luminosity drift away from the observations. 
Thus, we can also discard these models. 

Based on the above observations we tested other parameter combinations, as 
shown in Figure 15.161 In the upper panels we show the tests of models with slightly 
smaller energies and larger radii as compared with the reference model, and two 
different degrees of ^^Ni mixing. While the match of the LC for these models is 
satisfactory, the problem with the early-time velocities reappears, so we discard 
these models. In the lower panels of Figure 15.161 we show further combinations of 
parameters. Here we obtain an improvement in the LC while not compromising the 
agreement in the velocities. Among these models the one called ]Vi14r10e09ni56m 
provides the best match to the LC and velocity data. We therefore conclude that 
the 14 Mq scenario for SN 1999em cannot be ruled out, although the comparison 
with the data is not as good as that of the 19 Mq model (see Figure [573] and [5^ . 
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Figure 5.15: Comparison between models (lines) and observations (points) of 
SN 1999em for a low value of the pre-supernova mass of M = 14 Mq and the 
same mass of ^^Ni as for the reference model m19r8e1.25ni56. (Left panels): 
bolometric LCs. (Right panels): photospheric velocity evolution. (Upper pan- 
els): models m14r8e1ni56, m14r10e08ni56 and m14r12e08ni56 (see TablelO])- 
(Lower panels): models with different mixing of ^^Ni, E = 1 foe and R = 800 i?©. 
Note that the degree of ^^Ni is indicated in each panel as a fraction of the initial 
mass of the model. 
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Figure 5.16: Comparison between models (lines) and observations (points) of 
SN 1999em for a low value of the pre-supernova mass of M = 14 Mq and the 
same mass of ^^Ni as for the reference model m19r8e1.25ni56. (Left panels): 
bolometric LCs. (Right panels): photospheric velocity evolution. (Upper pan- 
els): models m14r8e1.1ni56I, m14r10e09ni56I and m14r10e09ni56M. (Lower 
panels): models m14r12e09ni56M and m14r12e09ni56I (see Table [53]). 
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Table 5.3: Model Parameters 
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The degree of ^^Ni mixing is given as a fraction of the initial mass of the model 
(Mo) 
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Chapter 6 



Observed and Physical Properties 



In this chapter I use the database presented in chapter [2] to examine the observed 
and physical properties of this sample of SNe II-P. I begin by calculating bolometric 
light curves and effective temperature evolution for the complete sample of SNe II- 
P in § 16.11 Then, I proceed to derive different parameters that characterize the 
bolometric LC, and I find that the SN sample shows a range of 1.15 dex in plateau 
luminosity, and plateau durations between 64 and 103 days. Comparing the shape 
of the transition between the plateau and the radioactive tail, I find that the size 
of the drop in luminosity is in the range 0.35 to 1.46 dex. The radioactive nickel 
masses calculated for all the SNe in our sample except for SN 1992am are below 
0.1 Mq. In section [672] I calculate a grid of hydrodynamical models for a range of 
initial masses, radii, explosion energies, nickel masses, and mixing. From this set 
of models and the observations I study correlations between different observed and 
physical parameters in § 16.31 
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6.1 Application to SN II-P Data 



In this section, we apply the cahbrations for bolometric corrections and effective 
temperatures presented in chapter § H] to our sample of 37 SNe II-P. Note that in 
order to calculate bolometric light curves we just need to apply equation fl4.4l) of 
§ I4.1.2[ Thus, we need to have (1) BVI photometry, (2) extinction corrections 
due to our own Galaxy, (3) host-galaxy ext inction correct i ons, a nd (4) distances. 
We obtained Galactic extinction values from Schlegel et al.l ( 1998 ) (see Table EH] of 



chapter [2]). Host- galaxy extinctions an d distances for most of the present sample 



were calculated by 



Olivares et al. 



(120091 ) and they are presented in Table [273] of • 



Figures I^TT] - 16 . 1 0] present the bolometric light curves and effective temperatures 
for the 37 SNe II-P in our sample. In order to place all the SNe in the same time 
scale and facilitate the comparison among them, the origin of time used is the middle 
point of the transition between the plateau and the radioactive tail [tpT] see § 12.31 
and Table [273]) . A simple inspection of the resulting bolometric hght curves reveals 
a high degree of heterogeneity and a variety of LC morphologies. 



There are ~16 SNe which were observed previous to the plateau phase, i.e. dur- 
ing the adiabatic cooling phase. Among these, there are three subluminous objects, 
SN 1999br, SN 2003bl and SN 2005cs which, in comparison with the rest, appear to 
have a steeper slope during the adiabatic cooling phase and a flatter plateau. Note, 
however, that in none of the cases we have data covering the maximum produced 
promptly after shock breakout, as predicted by the models. This lack of observation 
at maximum light is consistent with the prediction that such peak has a very short 
duration (a few hours), making it very difficult to observe. Only in recent years it 
was possible to observe this elusive phase and to detect the shock breakout for two 
SNe II-P thanks to observations at U V wavelength or to the grater cadence of t he 



the Supernova Legacy Survey (SNLS) (iGezari et al 



2008 



Tominaga et al 



200 



In all cases where we observe a plateau the bolomet r ic lum inosity remains nearly 
constant, slowly decreases, but never increases. lUtrobinI (]2007l ) has predicted, and we 
have also tested this with our models, that part of the dense core of the massive star 
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should be ejected in order to produce a plateau phase as observed. On the contrary, 
if the whole core remains in the compact remnant and the density profile encountered 
by the shock wave in the envelope is relatively flat, the resulting luminosity increases 
with time (see for example Figure 13161 of § I3.1.4p . Our data suggest that all SN II-P 
ejecta involve part of the dense core of the progenitor star. Also note that our sample 
of SNe shows different slopes — though never positive — during the "plateau" . Some 
SNe are more linearly declining instead of showing an obvious plateau in bolometric 
luminosity. Such is the case of SN 1999ca, SN 2002hx, SN 2002hj, SN 2003cn, 
SN 2003ci and SN 2003ip. We call these "intermediate-plateau" SNe. If we consider 
the definition of Type II-L SNe (SNe II-L) as those objects which show a linear, 
unint errupted luminosity decline until they reach the radioactive tail phase (iTuratto . 



20031 ) . then SN 2002hx could perhaps be classified as SN II-L. Interestingly, the 



existence of this group of intermediate SNe with tilted "plateaus" suggests a possible 
continuous connection between SNe II-P and SNe II-L, with the differences probably 
related to the progenitor mass. In this scenario, SNe II-L would have lower envelope 
masses probably due to mass lo ss during the proge nitor evolution, although there 



are alternative explanations (see ISwartz et al.l ( 1l99ll )). A quantitative classification 



based on th e average decline rate of the first 100 days of the light curve was given by 



Pat at et al. 



(jl994j ). In the following, we decided to include these intermediate SNe 
with the aim of showing how their observable properties compare with the genuine 
SNe II-P and we leave for future research a more precise definition of the limiting 
value of the "plateau" slope which would serve to divide both subclasses on the basis 
of their physical properties. It is important to keep in mind, however, that our code 
may not be suitable to model SNe II-L if such objects experience interaction with 
the circumstellar medium. 



To make a quantitative comparison between different SNe in our sample, we 
define a series of parameters that can be measured on the bolometric LC: (1) plateau 
luminosity {Lp), defined as the mean value between days —20 and —80 with respect 
to tpT', (2) plateau duration (Atp), defined as the interval of time (in the rest frame of 
the SN) during which the luminosity remains within ±0.2 dex of the mean value (Lp); 
(3) post-plateau drop (AlogL), defined as the logarithmic difference in luminosity 
between Lp and the tail luminosity calculated at day 20 with respect to tpx] and (4) 
^^Ni mass produced during the explosion (Mni). As examples we show in Figure 16.111 
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Figure 6.1: Bolometric luminosity (left) and effective temperatures (right) for 
SNe II-P as a function of time (part 1). The origin of tim e is the middle point 
between the plateau and the linear tail (tpT), as defined by lOlivares et al.l (120091 ) 
(see Table O of §[23D. 
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Figure 6.2: Bolometric luminosity (left) and effective temperatures (right) for 
SNe II-P as a function of time (part 2). The origin of time is the rn i ddle p oint 
between the plateau and the linear tail (tpr) as defined by Olivares et al.l ( 2009 ) (see 
Table Oof 
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Figure 6.3: Bolometric luminosity (left) and effective temperatures (right) for 
SNe II-P as a function of time (part 3). The origin of time is the rn i ddle p oint 
between the plateau and the linear tail (tpr) as defined by Olivares et al.l ( 2009 ) (see 
Table Oof 
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Figure 6.4: Bolometric luminosity (left) and effective temperatures (right) for 
SNe II-P as a function of time (part 4). The origin of time is the rn i ddle p oint 
between the plateau and the linear tail (tpr) as defined by Olivares et al.l ( 2009 ) (see 
Table O of §[23]) • 
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Figure 6.5: Bolometric luminosity (left) and effective temperatures (right) for 
SNe II-P as a function of time (part 5). The origin of time is the rn i ddle p oint 
between the plateau and the linear tail (tpr) as defined by Olivares et al.l ( 2009 ) (see 
Table Oof 
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Figure 6.6: Bolometric luminosity (left) and effective temperatures (right) for 
SNe II-P as a function of time (part 6). The origin of time is the rn i ddle p oint 
between the plateau and the linear tail (tpr) as defined by Olivares et al.l ( 2009 ) (see 
Table Oof 
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Figure 6.7: Bolometric luminosity (left) and effective temperatures (right) for 
SNe II-P as a function of time (part 7). The origin of time is the rn i ddle p oint 
between the plateau and the linear tail (tpr) as defined by Olivares et al.l ( 2009 ) (see 
Table Oof 
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Figure 6.8: Bolometric luminosity (left) and effective temperatures (right) for 
SNe II-P as a function of time (part 8). The origin of time is the rn i ddle p oint 
between the plateau and the linear tail (tpr) as defined by Olivares et al.l ( 2009 ) (see 
Table Oof 
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Figure 6.9: Bolometric luminosity (left) and effective temperatures (right) for 
SNe II-P as a function of time (part 9). The origin of time is the rn i ddle p oint 
between the plateau and the linear tail (tpr) as defined by Olivares et al.l ( 2009 ) (see 
Table Oof 
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Figure 6.10: Bolometric luminosity (left) and effective temperatures (right) for 
SNe 11-P as a function of time (part 10). The origin of time is the ni i ddle p oint 
between the plateau and the linear tail (tpr) as defined by Olivares et al.l ( 2009 ) (see 
Table O of §|23D. 



the bolometric LC for SN 1992ba, SN 1999em, SN 2002gw and SN 2003hl together 
with thier parameters (l)-(4). In Table [6?T] we list the values of the parameters and 
their uncertainties for for all the SNe in our sample. In the following subsection 
we give a detailed explanation of the calculation of these parameters, and their 
characteristic values and distributions. 



104 




40.5 



-100 -50 50 

days since transition time (tp.^) 



43 



42.5 



M 42 
J 41.5 

O 
J 

41 



40.5 



1 1 1 1 1 1 1 1 1 1 1 


SN 


1 1 1 1 _ 
1999em I 




















• 








• 








AL 


~ < < 1 , , , , 1 < < < 




, , 1 , " 



-100 -50 50 

days since transition time (tp.^) 



43 



42.5 



Dj 42 
J 41.5 

O 
J 

41 



40.5 



1 1 


1 1 1 1 


1 1 1 1 1 1 


SN 


1 1 1 1 _ 
2002gw ; 


















Lp 








I — 1 

• 














AL Z 


1 1 


1 , , , 


1 1 1 1 1 


1 1 


, , 1 , " 




-100 -50 50 

days since transition time (tp^) 



40.5 



-100 -50 50 

days since transition time (tp^) 



Figure 6.11: Bolometric LC and observed parameters [Lp, Atp, AL and Mni) for 
four SNe in our sample (see Table 16. ip . 
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Table 6.1: Observed Parameters for SNe II-P 





'o 
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LogLp 
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^ We generally use to(-E96) (see Table [231 of ij 12.31 1 as the explosion time. In case that this value is not 
available we use to(D05) (see Table [231 of S 12. 3p . If neither value is available we adopt estimations of to 
found in the literature (see 12.31 1. Note that the explosion time is used to calculate A/nj and both to{E96) 
and to{D05) provide consistent values for this parameter. 

*" Lower bound for Atp. 
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6.1.1 Plateau Luminosity: L 



In order to estimate the value of this parameter we average the bolometric luminosity 
between two well-defined epochs. The choice of —20 and —80 days for this was based 
on the behavior shown by most SNe. Column 4 of Table [HIT] shows the resulting values 
of this parameter and their uncertainties. The estimation of the uncertainty was done 
considering the one in the bolometric luminosity whose sources are (1) photometry, 
(2) extinction (see Table 1231) . (3) distance (see Table l2T3ll and (4) calibration of the 
bolometric correction. For the latter component, a value of 0.11 mag was adopted 
during the plateau, and of 0.02 mag during the tail (see § I4.1.2p . Figure I6TT2] shows 
the distribution of this parameter among our SNe. We see that the intermediate- 
plateau SNe have plateau luminosities which are consistent with the genuine SNe II-P. 
The histogram reveals that most of the SNe II-P (22 out of 37) have a characteristic 
plateau luminosity between 7.9 x 10^^ erg s~^ and 1.6 x 10^^ erg s~^. There is only one 
object, namely SN 1992am, with Lp > 3.16 xlO^^ g^g s"^ and two SNe (SN 1999br 
and SN 2003bl) with Lp < 5 xlO^^ erg s'\ 



1 ' ' ' ' I ' ' ' ' I ' ' ' 

No. of SNe: 37 




Figure 6.12: The distribution of Lp for SNe II-P. The SNe considered as intermediate 
between plateau and linear (see § 16.11) are marked in red. 
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In Figure EUS] we compare six exemplary SNe to illustrate that S Ne II-P display 
a wide range of plateau luminosities similar to the result reported by (jHamuyl . l200ll ) 
using V-band light curve. The range of luminosities encompassed by our sample is 
1.15 dex, equivalent to more than one order of magnitude of spread in their radiative 
energy output. The weighted average value of the characteristic plateau luminosity 
in our sample is < Lp >= 1.26 ± 0.019 x 10^^ erg s~^. 



43 



42.5 



CD 



o 



42 



41.5 



o 
o 
o 




-1 I I r 

O SN 1992am 

® SN2003hg 

® SN 1999em 

O SN 1992ba 

® SN 2005CS 

O SN 1999br 



O o 



100 -50 

days since transition 



Figure 6.13: Bolometric light curves for five SNe of our sample showing the range of 
variation of plateau luminosities. 
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6.1.2 Plateau duration: At 



The estimate of this parameter requires the calculation of two characteristic times: 
(1) the point where the L — Lp + SL called t+, and (2) the point where the 

L — Lp — 5L (I/-), called Here, SL represents the range of luminosity adopted 
to define the plateau phase. We chose a value of 0.2 dex which is equivalent to 
a variation of 1 mag. Note that t+ always occurs before t_ given the monotonic 
behavior of the LC. Thus by definition, Atp — {t- — t+)/{l + z), where z is the 
redshift of the SN. Although it is straight- forward in theory to define Atp , in 
practice its calculation involves some problems especially for those SNe which do not 
have good temporal coverage. We found the following cases: (a) SNe with data for 
L > and L < (b) SNe with data only for L < L^, and (c) SNe with data 
only for L > L^. 

In most of cases (a), we can determine t+ and t- and their uncertainties by 
interpolating between the nearest points on both sides of and L_, respectively. 
When such points are separated by more than 30 days from each other an extrap- 
olation is performed based on the two points right below L+ for t+, or above L_ 
for In cases (b) and (c) it was only possible to find t+ or respectively, by 
extrapolation. 

In order to compute Atp for SNe with poorer coverage, we use data of SNe with 
the best coverage to derive typical time intervals which in turn allow to obtain i_ or 
t+. For the estimation of t-, we use the typical interval found between this moment 
and the transition time, i.e. At- — {tpT — t-)/{l+z). A weighted average, < At- >, 
and its uncertainty is computed for the 13 SNe which allow a direct determination 
of t-. We find < At- >— 8.3 ± 0.5 days, i.e, 8.3 days before the middle point of the 
transition between plateau and tail. This value is then used to obtain t- for other 
SNe based on their values of tpT- 

A similar procedure can be followed for the estimate of t+. In this case its 
apphcation requires knowledge of the explosion time (to). For 13 SNe with both 
estimates of t+ and to, we compute the interval between those times, At+ — (i+ — 
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to)/(l + z), and its uncertainty. A weighted average of these values is obtained, 
< At+ >= 22.0 ± 1.5 days. This means that, on average, SNe II-P enter the plateau 
phase approximately 22 days after explosion. This result is in very good agreement 
with predictions from hydrodynamical models. We are thus able to estimate t+ 
using the average above only when the explosion time is known. There are ten 
additional objects for which we can only estimate a lower limit of t+ and thus of 
Atp. Column 3 of Table 16.11 shows the values of Atp and their uncertainties. The 
values corresponding to lower bounds of Atp are preceded by the > symbol. 

Figure [6.141 shows the distribution of Atp (excluding lower limits). The his- 
togram shows that most SNe (20 out of 27) have Atp values between 75 and 105 
days and the peak is around 90 days. Although the sample is small, we see hints of 
a bi-modal distribution with a secondary peak lying at about 60 days. However note 
that all but two of the SNe that form this secondary peak are intermediate-plateau 
SNe. 




40 60 80 100 

Atp [days] 



Figure 6.14: Distribution of plateau durations, Atp, values for SNe II-P. The SNe 
considered as intermediate between plateau and linear (see § 16. ip are marked in red. 
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Figure 16.151 shows a comparison of the hght curves of five SNe selected to 
illustrate the variety of plateau durations. The range of plateau durations is between 
47 and 103 days. However, if we restrict the sample to those SNe with bolometric 
LC showing a genuine plateau, the range is reduced to 64 - 103 days. The weighted 
average value of the plateau duration in our whole sample is < Atp >= 90.47 ± 1.18 
days. 
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Figure 6.15: Bolometric light curves for four SNe in our sample showing the range 
of variation of the plateau lengths. 
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6.1.3 Luminosity drop: AlogL 



The purpose of this parameter is to quantify the luminosity drop between the plateau 
and the radioactive phase which is sensitive to the nickel mass and its level of mixing. 
The calculation of AlogL involves the use of the plateau luminosity Lp (see I6.1.1"]) . 
and an estimation of the tail luminosity (L^) at some specific moment during the 
radioactive tail phase. The luminosity drop is simply AlogL = log Lp — logLj. 



We choose day 20 with respect to tpx (jOlivares et al.L l2009l ) as a convenient 
moment for the calculation because at that time all SNe have fully in the radioactive 
tail phase. It is therefore necessary to have at least one data point during the tail 
phase. In cases where more than two data points are available, a straight-line fit 
to the data is used to provide the value of Lt and its uncertainty. If only two data 
points are available, an interpolation is used. We were able to perform such fits and 
interpolations for 20 SNe, which additionally allowed us to determine a characteristic 
slope during the tail phase, < slope >= —0.004 dex day^^ which is in agreement 
with the slope expected for radioactive decay of ^^Co. This slope is used to derive 
Lt for those SNe which only have one data point during the radioactive tail. 



The values and uncertainties of Lj are given in column 5 of Table 16.11 Fig- 
ure [6]T6] shows the distribution of AlogL with a sharp peak at AlogL = 0.75 and 
a concentration of SNe (18 out of 26) in the range between 0.7 and 1. Note that the 
intermediate-plateau SNe have A log L values consistent with the genuine SNe II-P, 
with the exception of SN 2003cn that shows the largest value of A log L as a con- 
sequence of its low estimated ^^Ni mass (see next section). Also note that there is 
only one object, namely SN 2002fa, with AlogL < 0.4. 
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Figure 6.16: Distribution of A log L for SNe II- P. The SNe considered as intermediate 
between plateau and linear (see § 16.11) are marked in red. 
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In Figure EiTT] we compare the shape of the transition between the plateau phase 
and the radioactive tail for SNe having a well-sampled transition. The behavior of 
the transition is related to the ^^Ni mass synthesized in the explosion and its degree 
of mixing. M ore ^^Ni mixing produces a mo r e grad ual transition between the plateau 



and the tail ( lEastman et al 



1994; 



Utrobin, 



20071). The drop for SN 2003gd appears 



to be steeper and larger than for the other SNe, which is indicative of less mixing 
and a smaller ^^Ni mass. For the whole sample the luminosity drop is in the range of 
0.35-1.46. The weighted average value of the luminosity drop for the current sample 
is < A log L >= 0.783 ± 0.054. 
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Figure 6.17: Bolometric light curves of six SNe of the sample with different transition 
properties. 



114 



6.1.4 ^^Ni mass: Mm 

The bolometric luminosity during the radioactive tail phase is a direct indicator of 
the mass of ^^Ni synthesized during the explosion if we assume that all of the 7-rays 
resulting from ^^Co — > ^^Fe are fully thermalizecjl]. This a reasonable assumption 
given that most of the SNe II-P have a late-time decline rate consistent with ^^Co — > 
^^Fe decays (see § I6.1.3p . at least in the initial phase of the nebular era. However, if 
some 7-rays escape before being thermalized, the estimate of M-^i must be considered 
as a lower limit. To calculate MNi, we assume that gamma-ray deposition (see § 13.1.3^ 
is equal to 1 in the whole structure of the object and thus Lhoiit) = ei.ad(^) during 
the tail, with erad(i^) given by equation fl3.12p of § 13.1.31 Therefore, Mni in units of 
Mq is given by 



\M.] = Lhoijt) 

®^ 6.41 X 10^3 exp(^)+ 1.33 X 10^3 exp(^^)' 



(6.1) 



where TNi = 8.8 days, and tco = 113.6 days are the mean lifetimes of the radioactive 
isotopes. Note that to is the explosion time. Based on equation (16. ip we can calculate 
for the SNe with known to and tail luminosity. Specifically, we calculate a value 
of and its uncertainty for each data point during the tail. The resulting M^i is 
the weighted average of these values. Given that for some SNe we have two estimates 
of to ito{E96) and to{D05) see § 12. 3p . two different values of Mni can be calculated. 
However, in column 6 of Table 16.11 we present only one value of Mni. This is in 
general the value corresponding to to{E96). In cases that this value is not known, 
we adopted to{D05). And if neither value is known we used estimates of to found in 
the literature (see § l2.3p . We remark that the values of M^i calculated using to{D05) 
are systematically larger but consistent within the uncertainties with those based on 
to{E96). That is why we present only one value of M^i. 

Additionally, we include in our analysis SNe with t+ (see § I6.1.2p well deter- 
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Co is daughter of ^^Ni, which has a half-hfe of only 6.1 days 
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mined. Using and the characteristic < At^ > (see § I6.1.2p we can estimate the 
explosion time, t^, and therefore calculate Mni for more SNe in our sample. The un- 
certainty associated with this calculation is large because M^i is strongly dependent 
on the explosion time and the calculation of the bolometric luminosity (see § I6.1.ip . 

Figure 16.181 shows the distribution of Mni for all the SNe for which we can 
estimate this parameter. We do not include in this analysis SN 1992am for which 
we find a significantly higher value of the ^^Ni mass as compared with the rest of the 
sample (~ 0AM q see Table [6?T|) . Considering that there is a wide range of values 
for the explosion time of this SN in the literature, we checked if this high value of 
Mni could be due to an uncertainty in the value of to assumed. For example, if we 
use the estimate of to derived from t+ and < At+ > we find that the Mni is reduced 
to 0.26 ± 0.07 Mq. Although smaller, this value is still significantly large, which 
may imply that there are SNe II-P than can produce significant nucleosynthesis of 
radioactive material. Nevertheless we decided to exclude this SNe in the following 
analysis and leave it for a closer scrutiny in the future. The histogram shows that the 
SNe in our sample produce less than O.IMq of ^^Ni, with a peak at 0.035Mq. There 
is a slight tendency to a bi-modal distribution, but this should be confirmed with a 
larger sample of SNe. With the exception of SN 2003cn that shows the lowest value of 

in the whole sample, the intermediate-plateau SNe have values of < O.IM© 
consistent with the genuine SNe II-P. 

In Figure 16.191 we compare the tail luminosity for a subset of five SNe. Note that 
the amount of ^^Ni synthesized by the SN determines the height of the radioactive 
tail. The weighted average of ^^Ni masses obtained in our sample of SNe II-P is 
Mn; = 0.024 ±O.OO4M0. 

Before proceeding to the next section, we note that with the exception of 
SN 2003cn we did not find any peculiarity in the values of the Lp, AlogL and 
Mni for the intermediate-plateau SNe included in the analysis. This consistency 
strengthens the idea of a connection between linear and plateau subtypes. 
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No. of SNe: 19 




mediate between plateau and linear (see § I6.ip are marked in red. 



117 



t3f) 
0) 



o 



T 1 



o 
o 



o 



m 



o 



QD 



%1 



o 



Oo 



# 



-1 I I I l— 

O SN 1992af 
« SN 1999em 
m SN 2003hn 
® SN 2004dj 
• SN 2003gd 



O 



m 



d 



50 100 

days since transition 



Figure 6.19: Bolometric light curves of five SNe of our sample with different tail 
luminosities. 
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6.2 Grid of Hydro dynamical Models 



With the aim of comparing parameters between models and observations, and of 
studying correlations between observables and physical parameters, we calculated a 
grid of hydrodynamical models for different values of the initial mass (Mq), initial 
radius {Rq), explosion energy (E), ^^Ni mass {M^i) and mixing of nickel. Other 
parameters like chemical composition and density profile shape were kept fixed. The 
ranges of physical parameters used were: Mq = 10, 15, 20 and 25 Mq, E = 0.5, 1, 2 
and 3 foe, Rq = 500, 1000, and 1500 R©, Mm = 0.02, 0.04 and 0.07 Mq, and mixing 
of nickel out to fraction of Mq of 0.23, 0.5 and 0.8. 

For each model we calculated Lp, Atp, and AlogL using the same definitions 
as in the previous section to perform a consistent comparison between model and ob- 
servation. Additionally, we calculated the photospheric velocity at —30 days (f_3o) 
with respect to tpT- In Table [62] we present our results for a total of 46 hydrody- 
namical models. Although the table includes all the models, some sets of physical 
parameters seem unlikely to happen in nature. For example, models with initial 
masses of 15, 20 and 25 Mq and explosion energy of 0.5 foe, independently of the 
values of Rq, Mni and nickel mixing, yield values of Atp > 110 days, which is outside 
the observed range for our sample of SNe. Thus, we can rule out models and refine 
the set of physical parameters that could produce SNe II-P. 



In Figure 16.201 we show the sensitivity of our model light curves to the varia- 
tion of physical parameters for some selected models of Table 16. 2[ A glance at this 
figure indicates that (a) increasing the injected energy results in higher luminosi- 
ties and shorter plateau durations, (b) more massive progenitors produce longer 
plateaus and lower plateau luminosities, (c) more extended progenitors produce 
higher plateau luminosities and longer plateaus, and (d) a greater amount of ^^Ni 
produces a longer plateau and a higher tail luminosity. The behavior of our bolomet- 



ric LC on variation of physica. 



works (ILitvinova fc Nadezhin 



parameters i s consis t ent w ith those shown in previous 



1983 



1985 



Popovl . 1 19931 . among others). Based on 



the results of our models we study in greater detail the dependence of observable 
parameters on physical quantities {E, Rq, Mq and M^i). Figures [6.211 16.221 16.231 
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t [days] t [days] 

Figure 6.20: Bolometric LC for different (upper left) injected energies (E), (lower 
left) initial radius (Rq), (upper right) initial mass (Mq), and (lower right) ^^Ni 
mass (Mni). The shape of the LC, being sensitive to the hydrodynamics, is a useful 
tool to infer physical parameter of the SN progenitor. 

and 16.241 show Lp, Atp, AlogL and f(— 30), respectively, versus explosion energy. 
In these figures, the size of the symbols is proportional to the progenitor mass and 
the shape is related to the initial radius (square symbols for _Ro = SOOi?©, circles for 
i?o = 1000i?Q and triangles for Rq = ISOOi?©). The colors indicate different values 
of ^^Ni mass (red color for = O.O2M0, blue for = O.O4M0 and cyan for 
= O.O7M0). In all these figures we only show models with a mixing of ^^Ni of 
0.5Mo. 

From Figure 16.211 it is clear that there is a strong correlation between plateau 
luminosity and explosion energy. Larger explosion energies produce higher plateau 
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Table 6.2: Model Parameters 



Model 


Energy 
[foe] 


Mass 

[Mq] 


Radius 

[^0] 


If 


mixing'^ 


Log Lp 
[erg s"^] 


Aip 

[days] 


ALogL 
[erg s"^] 


„b 
[km s-l] 


1 


0.50 


15 


1000 


0.04 


0.23 


41.70 


119.91 


0.519 


1213 


2 


1.00 


15 


1000 


0.04 


0.23 


42.02 


86.16 


0.711 


2093 


3 


2.00 


15 


1000 


0.04 


0.23 


42.39 


62.48 


0.992 


3531 


4 


3.00 


15 


1000 


0.04 


0.23 


42.61 


45.71 


1.167 


4809 


5 


0.50 


15 


1000 


0.02 


0.5 


41.72 


97.79 


0.773 


1534 


6 


1.00 


15 


1000 


0.02 


0.5 


42.07 


78.80 


1.016 


2479 


7 


2.00 


15 


1000 


0.02 


0.5 


42.44 


53.02 


1.310 


4023 


8 


3.00 


15 


1000 


0.02 


0.5 


42.63 


45.07 


1.466 


5226 


9 


0.50 


15 


1000 


0.04 


0.5 


41.76 


116.43 


0.536 


1471 


10 


1.00 


15 


1000 


0.04 


0.5 


42.07 


87.83 


0.717 


2369 


11 


2.00 


15 


1000 


0.04 


0.5 


42.43 


64.13 


0.995 


3815 


12 


3.00 


15 


1000 


0.04 


0.5 


42.63 


48.07 


1.155 


5080 


13 


0.50 


15 


1000 


0.07 


0.5 


41.86 


135.94 


0.426 


1413 


14 


1.00 


15 


1000 


0.07 


0.5 


42.08 


98.59 


0.507 


2276 


15 


2.00 


15 


1000 


0.07 


0.5 


42.41 


74.19 


0.739 


3628 


16 


3.00 


15 


1000 


0.07 


0.5 


42.62 


59.03 


0.912 


4943 


17 


0.50 


15 


1000 


0.04 


0.8 


41.82 


112.78 


0.560 


1565 


18 


1.00 


15 


1000 


0.04 


0.8 


42.11 


84.85 


0.738 


2543 


19 


2.00 


15 


1000 


0.04 


0.8 


42.46 


60.91 


1.013 


4026 


20 


3.00 


15 


1000 


0.04 


0.8 


42.64 


49.89 


1.168 


5384 


21 


0.50 


20 


1000 


0.04 


0.5 


41.73 


151.32 


0.594 


1244 


22 


0.50 


25 


1000 


0.04 


0.5 


41.70 


177.66 


0.693 


1093 


23 


1.00 


10 


1000 


0.04 


0.5 


42.23 


43.13 


0.747 


3170 


24 


1.00 


20 


1000 


0.04 


0.5 


41.98 


116.04 


0.703 


1930 


25 


1.00 


25 


1000 


0.04 


0.5 


41.94 


138.93 


0.768 


1660 


26 


2.00 


20 


1000 


0.04 


0.5 


42.29 


92.64 


0.922 


2985 


27 


2.00 


25 


1000 


0.04 


0.5 


42.22 


110.99 


0.950 


2519 


28 


3.00 


20 


1000 


0.04 


0.5 


42.47 


82.43 


1.084 


3860 


29 


3.00 


25 


1000 


0.02 


0.5 


42.39 


99.09 


1.124 


3260 


30 


1.00 


10 


1000 


0.02 


0.5 


42.24 


37.21 


1.100 


3519 


31 


1.00 


20 


1000 


0.02 


0.5 


41.96 


108.36 


0.982 


1923 


32 


1.00 


25 


1000 


0.02 


0.5 


41.90 


131.66 


1.143 


1655 


33 


1.00 


10 


1000 


0.07 


0.5 


42.22 


65.77 


0.546 


2889 


34 


1.00 


20 


1000 


0.07 


0.5 


42.03 


125.44 


0.526 


1887 


35 


1.00 


25 


1000 


0.07 


0.5 


41.10 


147.83 


0.585 


1668 


36 


0.50 


15 


500 


0.04 


0.5 


41.64 


126.58 


0.402 


1373 


37 


1.00 


15 


500 


0.04 


0.5 


41.87 


92.70 


0.498 


2280 


38 


1.00 


15 


1500 


0.04 


0.5 


42.18 


83.05 


0.855 


2404 


39 


2.00 


15 


500 


0.04 


0.5 


42.23 


70.20 


0.770 


3656 


40 


2.00 


15 


1500 


0.04 


0.5 


42.54 


57.53 


1.130 


3877 


41 


3.00 


15 


500 


0.04 


0.5 


42.41 


59.36 


0.933 


4997 


42 


3.00 


15 


1500 


0.04 


0.5 


42.73 


47.90 


1.279 


5180 


43 


1.00 


15 


500 


0.02 


0.5 


41.86 


82.41 


0.766 


2381 


44 


1.00 


15 


1500 


0.02 


0.5 


42.19 


73.60 


1.163 


2516 


45 


1.00 


15 


500 


0.07 


0.5 


41.92 


106.16 


0.338 


2124 


46 


1.00 


15 


1500 


0.07 


0.5 


42.18 


93.49 


0.623 


2309 



^ The degree of ^^Ni mixing is given as a fraction of the initial mass of the model (Mo) 

is the photospheric velocity at —30 days with respect to the transition between the plateau and radioactive 
tail phases. 



121 



luminosities (see also upper-right panel of Figure [^■20p . However, there is significant 
dispersion in this relation (~0.4 dex). Note that for a fixed value of explosion en- 
ergy, more massive progenitors produce lower plateau luminosities and smaller initial 
radius also produce lower values of Lp. There is an additional contribution to the 
dispersion due to the ^^Ni mass, although it is lower than for the other parameters. 




Figure 6.21: The characteristic plateau luminosity {Lp) for all of the models pre- 
sented in Table 16.21 (with ^^Ni mixing of 0.5 Mq) as a function of explosion energy. 
The size of the symbols is proportional to the initial mass. The shape is related to 
the initial radius (square symbols for Rq = 5OOi?0, circles for Rq = lOOOi?©, and 
triangles for Rq = ISOOi?©). The colors indicate different values of ^^Ni mass (red 
color for Mm = 0.02Mq, blue for Mm = O.O4M0, and cyan for Mm^ = O.OTMo). 



The plateau duration also correlates with explosion energy (see Figure I6.22p 
but in this case the dependence is wearer than in the case of the plateau luminosity. 
The progenitor mass appears to be the most important factor affecting Atp. Note, 
for example, that almost the full range of variation of Atp is covered for fixed E = 1 
foe, only by changing the initial mass between 10-25 Mq. Neverthe less the effect of 
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^^Ni mass is not negligible. Increasing the ^^Ni mass the plateau duration becomes 
noticeably larger. The effect of ^^Ni mass seems to depend on energy, being greater 
for smaller values of explosion energy. The initial radius produces only a minor 
effect on At p. Another point is that we found several models with Atp < 75 days, 
interestingly these models tend to a show a linear behavior consistently with the 
observations. 
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Figure 6.22: The plateau duration (Atp) for all of the models presented in Table [62] 
(with ^®Ni mixing of 0.5 Mq) as a function of explosion energy. The size of the symbols 
is proportional to the initial mass. The shape is related to the initial radius (square 
symbols for Rq = 500Rq, circles for Rq = lOOOi?©, and triangles for Rq = 1500/?©). 
The colors indicate different values of ^^Ni mass (red color for M^i = O.O2M0, blue 
for Mm = O.O4M0, and cyan for Mm = 0.07Mq). 



From Figure [6. 231 we see that the post-plateau drop AL is slightly dependent on 
E. The ^^Ni mass appears as the main factor in the determination of this parameter, 
as expected from the large effect of Mni on the tail luminosity. Lower values of M^j 
lead to reduced brightness in the tail and therefore to higher values of AL, and vice 



123 



versa. Also note that AL is not very sensitive to the initial mass but there is some 
dependence on the initial radius. 












2 
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Figure 6.23: The drop in luminosity (AL) for all of the models presented in Table [62] 
(with ^^Ni mixing of 0.5 Mq) as a function of explosion energy. The size of the symbols 
is proportional to the initial mass. The shape is related to the initial radius (square 
symbols for Rq = 500Rq, circles for Rq = IOOO-R0, and triangles for Rq = 15OOi?0). 
The colors indicate different values of ^^Ni mass (red color for Mni = O.O2M0, blue 
for Mm = O.O4M0, and cyan for = O.OTM©). 



In Figure 16.241 we plot the photospheric velocity at —30 days with respect to 
tpT versus explosion energy. The selection of this particular epoch to evaluate the 
photospheric velocity was based on the work of lOlivares et al.l (120091 ) which calibrate 
the standard- candle method for estimating distances using expansion velocities at 
this epoch. We found a strong dependence of v_^q on the explosion energy, as 
expected from the fact that the explosion energy is almost completely converted into 
kinetic energy a few days after the explosion (see § IS.l.ip . Also note that the value 
of f _3o depends on the progenitor mass in the sense that more massive progenitors 
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produce lower velocities for a fixed energy. There is also some dependence on ^^Ni 
at least for lower values of the progenitor mass but we do not find a dependence on 
initial radius. 



The analysis presented above is a qualitative demonstration of the correlations 
among observable and physical parameters based on a limited amount of models. 
Further work is planned to enlarge the set of models in order to provide quantitative 
relations. 




Figure 6.24: The photospheric velocity at —30 days with respect to tpT (f (— 30) ) for 
all of the models presented in Table 16.21 (with mixing ^^Ni of 0.5 Mq) as a function 
of the explosion energy. The size of the symbols is proportional to the initial mass. 
The shape is related to the initial radius (square symbols for Rq = 500Rq, circles for 
Rq = IOOO-Rq, and triangles for _Ro = 15OOi?0). The colors indicate different values 
of ^^Ni mass (red color for = O.O2M0, blue for = O.O4M0, and cyan for 
Mni = O.O7M0). 
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6.3 Observed and Modeled Correlations 



In this section, we use the parameters Lp, Atp, ALogL, and Mni, calculated both 
from the observations of our SN II-P sample (see Table 16.11) and from our grid of 
hydrodynamical models (see Table \^72\i to analyze correlations among them. We also 
include in this analysis the expansion velocity r neasured from the Fe II A5169 line at 
—30 days with respect to tpT as calculated by lOlivares et al.l (120091 ) (see Table 4.1 
of that work). 

We begin by analyzing how is related to other observable parameters. 
Figures [6.251 16.26] 16.271 and 16.281 show Lp, Atp, ALogL and f-30, respectively, as a 
function of the ^^Ni mass. In these figures, filled symbols represent the models and 
open symbols show the observations. Different colors are used for identifying models 
with different explosion energies (magenta color for E = 0.5 foe, blue for E = 1 foe, 
green for = 2 foe, and cyan for E = 3 foe), and the size of the model symbols 
is proportional to the progenitor mass. We use . The intermediate-plateau SNe are 
shown with a different symbol (triangles) than genuine plateau SNe (squares). 



Hamuyi ( 120031 ) found a correlation between plateau luminosity and ^^Ni mass 
in the sense that more luminous supernovae produce more ^^Ni. The observational 
data in Figure [6.251 shows a slight trend in agree r nent w ith such correlation, although 



less pronounced than the one found by [Hamuyi ( l2003[ ). The models show a range of 



plateau luminosities for any specific value of ^^Ni mass which is mainly caused by 
the explosion energy. It is important to note that the ^^Ni mass is introduced as a 
free parameter in our model while it is actually determined by the nucleosynthesis 
produced during the explosion and by how much ^^Ni falls back on the compact 
remnant. In any case, we see that the observations are consistent with models of 
explosion energy ~1 foe and different progenitors properties (mass an radius). 

As shown in Figure I6.26[ the data show no correlation between Atp and ^®Ni 
mass. Thus is consistent with the models which show only a mild dependence of the 
plateau duration with ^^Ni mass. 
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Figure 16.271 shows that there is a strong correlation between A Log L and ^^Ni 
mass, both in the models and observations. The lower the of ^^Ni mass the greater 
the luminosity drop. 



Finally, Figure [6.281 compares ^^Ni mass with f_3o. The observations suggest a 
trend where fa s ter SN e produce larger amounts of ^^Ni, as was previously pointed 
out by iHamuyl (120031 ). Since f_3o is closely related with the kinetic energy which is 
a considerable percentage of the explosion energy, this correlation suggests that SNe 
with higher explosion energy produce more nuclear burning or suffer less fall-back of 
the reprocessed material. From the comparison between observations and models, 
we see that explosion energies around 0.5 foe occur, but are not realized very often 
in Nature. 



Given that the luminosity is an important parameter in the determination of 
distances, we examine the dependence of this parameter on other observables. In 
Figure I6.29[ I6.30[ and 16.311 we show its dependence on f_3o, Atp, and ALogL, 
respectively. The symbol shapes and colors have the same meaning as in the previous 
figures. 



It is known that there is a strong cor relation between 
pansion velocity during the plateau phase flHamuv fc Pintd. 



of SNe this correlation was studied in detail by lOlivares et al. 



umin osity and the ex- 



20021) . For our sample 



mm using -30 days 



with respect to tpT as the epoch where the correlation was calibrated. Figure 16.291 
shows this correlation for our models and the observational data. From the figure 
we see that the models reproduce very well the trend of the correlation, and we can 
identify the explosion energy as the main driving parameter. In addition, the mod- 
els show that the dispersion is related to the progenitor properties. For example, 
more massive progenitors produce lower expansion velocities for approximately the 
same luminosity. Although the trend of the correlation is well reproduced by the 
models we can see that there is a shift in the relationship between the models and 
the obse rvations. One r easoii for the shift could be related to luminosity calibration 
used by 



Olivares et al 



(|2009[ ) to estimate distances for the SNe in our sample. The 
zero point is based on two object, SN 1999em and SN 2004dj, with known Cepheids 
distance. These two SNe show a difference of ~0.9 mag in their SCM-calibrated 
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Figure 6.25: The characteristic plateau luminosity as a function of Ni mass for 
models (filled symbols) and observations (open symbols). The size of the symbols 
is proportional to the initial mass. Colors indicate the explosion energy: magenta 
for E = 0.5 foe, blue for E = 1 foe, green for E = 2 foe and cyan for E = 3 foe. 
The intermediate-plateau SNe are shown with triangles while squares are used for 
genuine plateau SNe. 
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Figure 6.26: The plateau duration as a function of ^^Ni mass for models (filled 
symbols) and observations (open symbols). The size of the symbols is proportional 
to the initial mass. Colors indicate the explosion energy: magenta for E = 0.5 foe, 
blue for E = 1 foe, green for E = 2 foe and cyan for E = 3 foe. The intermediate- 
plateau SNe are shown with triangles while squares are used for genuine plateau 
SNe. 
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Figure 6.27: The drop in luminosity as a function of ^^Ni mass for models (filled 
symbols) and observations (open symbols). The size of the symbols is proportional 
to the initial mass. Colors indicate the explosion energy: magenta for E = 0.5 foe, 
blue for E = 1 foe, green for E = 2 foe and cyan for E = 3 foe. The intermediate- 
plateau SNe are shown with triangles while squares are used for genuine plateau 
SNe. 
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Figure 6.28: The photospheric velocity at —30 days with respect to tpT as a function 
of ^^Ni mass for models (filled symbols) and observations (open symbols). The size of 
the symbols is proportional to the initial mass. Colors indicate the explosion energy: 
magenta for E = 0.5 foe, blue for E = 1 foe, green for E = 2 foe and cyan for E = 3 
foe. The intermediate-plateau SNe are shown with triangles while squares are used 
for genuine plateau SNe. 
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absolute magnitude, which could lead to a shift of ~0.2 dex in luminosity. Other 
reason for the discrepancy may be that we are using photospheric velocities from 
the models and line velocities of Fe II A5169 (f fcIi) from the observations. The right 
way to do the comparison would be to convert f peii into a photospheric veloc ity us- 
ing the calibration based on atmospheric models given by iJones et al.l (120091 ). The 
problem is that this calibration is only valid for v > 5000 km s~^ (see § 14. 2p and 
therefore is out of range for most of our data. Another possibility is that our set of 
models does not cover the appropriate range of physical parameters. In this sense, 
we note that models with lower values of the initial mass are more consistent with 
the observations. We will explore these issues in greater detail in a future research. 



From the models shown in Figure 16.301 we see that there is a relation between 
plateau duration and plateau luminosity. This relation is essentially a consequence of 
the explosion energy but we see an important dispersion associated with the progen- 
itor properties, especially the progenitor mass. The relationship is not present in the 
observational data although the uncertainties in the estimation of plateau durations 
are large and the statistics are small. We also see that models with the lowest values 
of explosion energy and hight masses are not consistent with the observations. 

Figure 16.311 shows that the values of Lp and AL from models and observations 
are in very good agreement, especially for models with E = 1 and 2 foe. We also see, 
for the models, a dependence between these two parameters which is not observed 
in the data. 



Finally, we compare the behavior of ti-so with respect to the plateau duration 
in Figure 16.321 We see that models predict a strong correlation between plateau 
durations and expansion velocities. This correlation can be understood based on 
the dependence of the expansion velocity on explosion energy. It is a well known 
fact that more energetic explosions produce shorter plateaus and vice versa. The 
observational data appear to follow this relation but the dispersion and the errors 
are large, making it difficult to draw a definitive conclusion. We note, however, that 
the observations rule out the lowest energy values. 

Based on the comparison of different observables between models and observa- 
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Figure 6.29: The characteristic plateau luminosity as a function of vso for mod- 
els (filled symbols) and observations (open symbols). The size of the symbols is 
proportional to the initial mass. Colors indicate the explosion energy: magenta for 
E = 0.5 foe, blue for E = 1 foe, green for E = 2 foe and cyan for E = 3 foe. 
The intermediate-plateau SNe are shown with triangles while squares are used for 
genuine plateau SNe. 
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Figure 6.30: The characteristic plateau luminosity as a function of plateau duration 
for models (filled symbols) and observations (open symbols). The size of the symbols 
is proportional to the initial mass. Colors indicate the explosion energy: magenta 
for E = 0.5 foe, blue for E = 1 foe, green for E = 2 foe and cyan for E = 3 foe. 
The intermediate-plateau SNe are shown with triangles while squares are used for 
genuine plateau SNe. 
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Figure 6.31: The characteristic plateau luminosity as a function of the drop in lu- 
minosity for models (filled symbols) and observations (open symbols). The size of 
the symbols is proportional to the initial mass. Colors indicate the explosion energy: 
magenta for E = 0.5 foe, blue for E = 1 foe, green for E = 2 foe and cyan for E = 3 
foe. The intermediate-plateau SNe are shown with triangles while squares are used 
for genuine plateau SNe. 
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tions we conclude in general models that with E = 0.5 foe do not occur very often 
discarded. Also models with E = 3 foe and M < 25Mq ot E = 1 foe and M > 20Mq 
are unlikely to explain the observations of our sample of SNe II-P. 
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Figure 6.32: The plateau duration as a function of f_3o for models (filled symbols) 
and observations (open symbols). The size of the symbols is proportional to the 
initial mass. Colors indicate the explosion energy: magenta for E = 0.5 foe, blue for 
E = 1 foe, green for E = 2 foe and cyan for E = 3 foe. The intermediate-plateau 
SNe are shown with triangles while squares are used for genuine plateau SNe. 
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Chapter 7 



Summary 



The main propouse of this thesis is the development of a one-dimensional, Lagrangian 
hydrodynamic code to model bolometric light curves of SNe II-P with the ultimate 
goal of analyzing the physical properties of an available large sample of SNe with 
highly-precise, well-sampled observations. In the first part of the thesis I focused 
on the development of the hydrodynamic code, the incorporation of the required 
micro-physics and the construction of the progenitors models. 

Since our code produces bolometric light curves and effective temperatures, it 
was also necessary to calculate such quantities from the observed photometry be- 
fore performing a comparison between models and observations. Hence, the second 
(Chapter Hj) of the thesis was devoted to deriving reliable calibrations for bolometric 
corrections (BC) and effective temperatures based on BVI photometry and appli- 
cable to SNe II-P. In this analysis I used data of three well-observed SNe 1987A, 
1999em, and 2003hn; and two series of atmosphere models by Eastman et al. (1996) 
and Dessart & Hillier (2005b). I found a tight correlation betewen BC and photo- 
metric colors which allows one to estimate bolometric luminosity with an uncertainty 
of 0.05 dex during the plateau phase. On the radioactive tail I found that the BC 
was independent of color with a value of —0.70 mag and a scatter of 0.02 mag based 
only on the behavior of SN 1987A. I noticed the importance of including L- and 
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M-band photometry in order to derive an accurate BC during this phase. 



Then, I studied in detail the performance of the code via its apphcation to one 
of the best observed SNe II-P, namely SN 1999em, obtaining a very good agreement 
with the observations using the following physical parameters E = 1.25 foe, M = 
19 Mq, R = 800 i?0, and = 0.056 Mq. In this analysis, I found that an extended 
mixing of ^^Ni is needed in order to reproduce a plateau as flat as that shown by 
the observations of SN 1999em. I note that the plateau phase, at least in the case of 
this SN, is powered by the energy deposited by the shock wave and released by the 
recombination process, plus some extra energy deposited by the radioactive material. 
I compa red our rusults \ yith t wo p revious hydrod ynamical studies of SN 1999em 
given by Baklanov et al. (2005) and UtrobinI ( 2007 ). Even if these studies are based 
on a more sofisticaded prescription of the radiative transfer, there is a very good 
agreement among the results obtained by the three works. In particular, the physical 
pa rameters obtained in o ur c alculations are intermediate between those estimated 
J2OO5I ) and lutrobini feoOTf ). 



by 



Baklanov et al. 



The remarkably good agreement of our model with the observations of SN 1999em, 
and the consistency found with previous hydrodynamical studies of this SN using 
more sophisticated code give us confidence in our attempt to model SNe II-P at least 
during the plateau phase and early evolution of the radioactive tail phase, in spite 
of the simplifications assumed in our calculation. 



The last part of this thesis was focused on the study of the intrinsic properties of 
our sample of SNe. First, I used the calibrations of BC and effective temperature to 
derive bolometric luminosities and effective temperatures for our complete sample of 
SNe II-P. In order to make a quantitative comparison among the objects, I calculated 
a set of parameters that characterize the bolometric LC: (1) plateau luminosity 
(Lp), (2) plateau duration (Atp), (3) the drop in luminosity between plateau and 
radioactive tail (AlogL) and (4) the nickel mass (Mni). 

I found that the SN sample shows a range of 1.15 dex in plateau luminosity 
with a weighted average value of < Lp >= 1.26 ± 0.019 x 10^^ erg s~^, and plateau 
durations between 64 and 103 days with a weighted average value of < Atp >= 
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90.47 ±1.18 days. Comparing the shape of the transition between the plateau and 
the radioactive tail, I found that the size of the drop in luminosity to be in the range 
of 0.35 to 1.46 dex with a weighted average of < A log L >= 0.783 ± 0.054 dex. The 
radioactive nickel masses calculated for all of the SNe in our sample are below O.IMq 
with a weighted average of 0.024 ± O.OO4M0. However, it is important to note that I 
have excluded the extreme case of SN 1992am from the analysis. For this SN I found 
Mni > 0.26 ± 0.07 Mq implying that there are SNe II-P which produce significant 
nucleosynthesis during the explosion. We leave a detailed analysis of this SN for 
a future work. SNe with a more linearly declining behavior during the "plateau" 
phase show consistent values of Lp, AlogL, and Mni, which strengthens the idea 
of a connection between linear and plateau SNe II. In addition, our models show a 
strong dependence of the Atp on the progenitor mass in the sense that progenitors 
with lower masses produce shorter Atp. This could mean that line ar SNe in gener al 
have less massive progenitors, as has been previously proposed (see iTurattd (l2003l )). 



Then I proceeded to calculate a grid of 46 hydrodynamical models for a range 
of initial masses (Mq), radii (-Ro), explosion energies (E), nickel masses (Mni), and 
nickel mixing. For each model I calculated Lp, Atp, A log L and the photospheric ve- 
locity (vph), and I studied the dependence of these observable parameters on physical 
quantities {E, Rq, Mq and M^i). 

From this analysis I found that the explosion energy is the main factor in the 
determination of Lp and Vph, with a smaller dependence on Mq, and on Rq for the 
case of Lp. I found that Mni has no noticeable effect on the those parameters. In 
the case of Atp, I identify the progenitor mass as the most important parameter but 
E and Mnj also produce an important effect. Finally, I found that AlogL largely 
depends on Mm, with some dependence on E and Rq but not on Mq. 

I studied correlations among the observable parameters, Lp, Atp, A log L, M^, 
and Vph calculated both from the observations and hydrodynamical models. The 
most important results of this analysis are: (1) our models reproduce the strong 
correlations shown by the observations between luminosity and expansion velocity 
during the plateau. The slope of the relation is consistent with the observations, 
but there is a general shift between observational and model data. The shift is 
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possibly related with the zero point of th e distance calibratio n used to calculate 
distances for the SNe in our sample (see lOlivares et al.l . 120091 ). which could lead 



to a shift of ~0.2 dex in luminosity. I identified the explosion energy as the main 
driving parameter of this correlation. (2) I did not find a clear correlation between 
Mni and other observables with the exception of A log L where I found that lower 
values of correspond to higher AlogL and vice versa. (3) Our models show a 
certain relationship between Lp and Atp in the sense that more luminous SNe have 
shorter Atp, but this is not seen in the observa tions. It is intresting to n ote that this 
correlation was also found in a recent work by iKasen &: Woosleyl ( 120091 ) using their 
hydrodynamical code and initial models from stellar-evolution calculations. (4) The 
models show a relationship between Atp and the expansion velocity with smaller 
values of Atp associated to SNe with greater expansion velocities. The observations 
show a similar trend but the scatter is large and the statistics is too poor to provide 
a definitive conclusion. 



Finally, based on the comparison of different observables between models and 
observational data we conclude in general models with (a) E = 0.5 foe and M > 
lOM©, (h) E = 3 foe and M < 25Mq, and (c) E = 1 foe and M > 20Mq do not 
occur very often. 
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Appendix A 

Finite Difference Approximations 



A.l Computational grid 

A finite difference representation of a system of time dependent equations, such as 
tlie one presented in § 13. 1.1 1 (equations 13. Ifl3.5p . consists of defining a space-time grid 
upon whicli finite difference approximations are performed to evaluate tlie functions 
and tfieir derivatives which constitute the partial differential equations under study. 
As an example, Figure lA.ll shows a schematic space-time grid which consists of a 
set of fixed coordinate positions Xk defined at discrete times t". Throughout this 
discussion, the subscript k represents the spatial index and n, the temporal index. 
The temporal variable changes in a discrete fashion such that t""*"^ = t"' + At"'"'"^/^, 
where At""^^/^ is the time step between t""^^ and t". If the spatial coordinates are 
fixed, as is the case of Lagrangian coordinates, they do not carry a temporal index. 
The interval between two spatial coordinates, Xk+i and Xk, is denoted by Axk+1/2 
and it is commonly referred to as grid cell or zone. Note that Axfc_(_i/2 and At"+^/^ 
need not be constant throughout the grid. 

Dependent variables are defined on this space-time grid and the evolution equa- 
tions take them from one spatial level to the next. However, there are different ways 
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to represent a dependent variable on the grid. For example, given the function u{x,t) 
this can be defined at cell corners denoted by = u{xk,t"'), at cell edges denoted 
by ulti^ = ^itj^ fn+i/2 = fn^ 0.5 Ar+^Z^ or at cell centers denoted 

by = M(xA;+i/2,r+^/^) with Xk+i/2 = Xk + 0.5 Axfc+i/2 (see Figure EI])- The 

choice of the representation is based on the physical meaning of each variable. 
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Figure A.l: Space-time grid. The values of the discretized function u{x,t) can be 
defined either at cell corners (e.g., u^), cell edges (e.g., ulZl^"^) or cell centers (e.g., 

n-l/2N 
^fc+1/2^ 



A. 2 The Difference Equations 



In the Lagrangian prescription of the hydrodynamic equations in spherical sym- 
metric approximation (equations I3.1f[33] of Section |3TT|) . the object is represented 
by concentric spherical shells (zones) of mass which represent the spatial 

coordinates of the system (analogous to Axfe+1/2 in the previous discussion). The 
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radii are defined at cell corners, the velocities v^'^^^'^ are defined at the cen- 
ters of vertical cell edges and the intensive quantities (p, P and E) are defined at 
the centers of the horizontal cell edges; for example the pressure is represented as 
^k+i/2- "^^^ time steps are also required: one to advance the velocity At" defined by 
At"= 0.5(At"~-^/^ + At""*"^/^), and another to advance the material state variables, 

The calculation procedure is to take known radii and velocities at and de- 
termine the new values f^^^^^, r^'^^ and VJ^^^^ values k using the following 
explicit difference representation of the hydrodynamics equations 



' k 



where Q is the artificial viscosity term included in the equation to take into account 



numerical discontinuities which are produced by the shock wave. For this 



;erm I 



have adopted the following expression given by IVon Neumann fc Richtmyerl (1l950l ) 



I ^ / n+1/2 _ n+l/2x2 n+1/2 n+1/2 n+1/2 

gn+1/2 ^ J '-'<3 y^k+l/2 ^k ) Pfc+1/2 1^ %+l ^ ^k 

k+1/2 ] g otherwise 



where Cg is a constant set to a value of 2 in our calculations. The time step At""*"^/^ is 
chosen subject to assure stability and accuracy. The stability of the explicit numerical 
scheme is obtained if the time step is constrained to the Courant condition, i.e.. 



143 



^ nim{Am^+}/2/<+i,/J, (A.4) 

where v'^'j^^-^^^ is the sound speed which varies from zone to zone. The time step 
is taken as a fraction of the Courant condition for all zones. Additional conditions 
on the time step are imposed to assure the accuracy of the results: it is required 
that changes in temperature, density and flux over one time step be less than 5%. 
This constraint added to that on the time step is a consequence of the iterative 
procedure used to solve the radiative transfer (see following discussion). Note that 
the iterations usually do not converge if At is too large. Even if we do not need to 
iterate, the size of At in the implicit scheme is limited by the error growth rate. To 
keep the error low, changes in the physical quantities over one time step must be less 
than, for example, 5%. So, the time step is limited, by reason of accuracy, to the 
physical time scale on which quantities change. However, in general, the physical 
time scale is much longer than the Courant time. 

The energy and radiation transfer equations (equation 13.41 and 13.51 of Sec- 
tion 13. ip are solved using an iterative procedure. The discretization of the energy 
equation in the diffusion limit for a semi-implicit scheme is written as 



~ ^k+1/2 + o(-Pfc+l/2 + -Pfe+1/2 + Qk+1/2 + Qk+1/2) ^Vk+1/2 



1 

k+l/2 ~ ^k+1/2 2'^-' fc+1/2 fe+1/2 ^k+1/2 ~^ ^k+l/2) 

Afn+1/2 

-{e{LiX\ - Lr^) + (1 - e){Li^, - m^;^^ (A.5) 



where AVfc+1/2 = '^k+1/2 ~ written by numerical convenience as 



/Am 



fc+1/2 
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Here 6' is a parameter with reflects the time-centering: 6 = 1 fully implicit scheme, 
and 6 = 0.5 for time-centering. This latter value for 6 is the one I usually adopted 
in the calculations. The luminosity is written as 



6K Am A; 



where = is the "fiux-limiter" , and is written as 



-Ki, 



Here, k, which also appears in equation \A.Q\ is the Rosseland mean opacity evaluated 
at rrik and t", i.e., = = k{T^,p^). Note, however, that the intensive quantities 
are known at the middle points of the grid. Therefore, it is needed to average the 
value of K between the points mk-1/2 and mk+i/2- Special care must be taken to do 
this average in order to prevent numerical noise which appears when the radiation 
flow is propagated through a steep front (su ch as a shock-wave front) where T and 



P change abruptly. Following IChristyl (119671 ). instead of /t, I have used an effective 
opacity defined by 



1 \ ^fc+l/2/^A:+l/2 + T^-i/2/'^k- 



1/2 



e// -^^+1/2 ^ -'fc-1/2 



Note that knowing the values of the temperature -^fc+1/2 -^fc+3/2 

values k, the transport problem is solved. However, equation (1A.5P is nonlinear 
with respect to the new temperatures at n + 1, according to the dependence of 
j^n+i^pn+i on these new values. Defining an extrapolated temperature 

as TJ^^y^ = T'^_^^i2 + <^^fc+i/2) one can linearize the energy equation in corrections 
5Tfc+i/2, which leads to a system of N (equal to the number of points grid) coupled 
linear equations. Each equation takes the form 
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Ak+1/2 STk+3/2 + Bk+1/2 5Tk+l/2 + Ck+1/2 = D 



k+1/2, 



(Ai 



where = G{Tk^i/2,Tk^i/2,Tk+i/2)\t=t", is a function of the temperature eval- 

uated at only three spatial points and at the previous time step t". The other 
constants in equation IA.8I are derivatives of this function "G" with respect to ST at 
the three different spatial points evaluated at t", i.e., 



A 



dG 



k+1/2 



di5Tk+3/2) 



B 



dG 



fe+1/2 



d{5Ty 



k+1/2} 



and G k+1/2 



dG 



1=1" 



diSTk^i/2) 



t=t^ 



Because each equation of the system depends on the temperature at only three 
adjacent mass points (see equation lA.Sp . the system of algebraic equations takes the 
from of a tridiagonal matrix. This tridiagonal system is solved for all the correction 
temperatures ST (inverting the tridiagonal matrix). The zone temperatures are 
corrected and the process is repeated until the corrections are reduced to a sufficiently 
small amount. The condition that is usually adopted in the code for the convergence 
of the iterative process is that max(^) < 10~^, where max(^) is the maximum 
value of the ^ in all the grid points. 
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Appendix B 

Solution of Gamma-ray transfer 



In this appendix I describe the solution of the 7-ray transfer in gray approximation 
assuming that the 7-ray opacity is purely absorptive and independent of the energy, 
and allowing a spherically symmetric distribution for the sources of the 7 photons — 
i.e., the radioactive material. From this s olution, the ene rgy deposited by 7-rays can 
be calculated becuase, by definition (see IClaytonl . the total energy absorbed 

by matter per unit of time and mass is 



dE 

- — = I IdQ, 
dm 



where the integral is done in all solid angles, / is the specific intensity and /t^ is 
the opacity of the 7-rays for which we adopted the value = 0.06 ?/e cm^ g^^ 



([Sutherland &: Wheeled . Il984l ) and ye is the number of electrons per baryon. Al- 
though our code solves the equations in spherical symmetry, in order to correctly 
take into account the energy deposited in a specific point, it is necessary to consider 
rays which come from different directions. Nevertheless, the axial symmetry is pre- 
served. Therefore, the previous integral in spherical coordinates and assuming axial 
symmetry can be written as 
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dE 

— = 2nK^ I I(r,e)sined9. (B.l) 
am 



mm 



As can be seen in the above expression, to calculate the energy deposited by 7- 
rays one should find the specific intensity in each point of the object by solving the 
transport of 7-rays. The radiative transfer equation in gray approximation is 



| = + (B^2) 



where / and 5* are the specific intensity and the source function integrated in fre- 
quency, and dr = p ds is the differential optical depth of 7-rays along the direction 
of the beam, s. The formal solution of this equation is 



J(r) = Jo e-" + f S{t) e'("-"')rfr'. (B.3) 
Jo 

The source function, S, is a known function in our problem and it depends on the 
distribution of the radioactive material. Specifically, it can be written as 



^ = e™d(t)e(r)K, (B.4) 

where eradif) denotes the time-dependent rate of energy released per gram by Ni- 
Co-Fe decay, 

e^ad = 3.9 X 10^° exp(-t/rNi) + 6.78 x 109[exp(-t/rco) - exp(-t/rNi)] erg g"^ s-\ 
and ^ (r) of equation IB. 41 is the initial distribution of ^^nI^. Note that this function 



^ Since ^^Co is a decay product of ^^Ni, its spatial distribution is described by the same function 
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